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LINEAR  ARRAYS:  CURRENTS,  IMPEDANCES,  AND  FIELDS,  III 

by 

Ronold  W.  P,  King  and  Sheldon  S.  Sandler 
Gordon  McKay  Laboratory  of  Applied  Science 
Harvard  University 
Cambridge.  Massachusetts 

ABSTRACT 

Further  results  are  given  for  the  curtain  array  with  three-quarter  and 
one-half  wavelength  elements.  Since  the  general  solution  assumes  an  indeter¬ 
minate  form,  a  special  solution  is  given  for  arrays  with  half -wavelength  ele- 
mients.  This  special  solution  is  obtained  from  the  original  Hallen  integral 
equation.  Calculations  are  given  for  two  three-element  endfire  arrays,  for 
three-quarter  and  one-half  wavelength  elements  and  the  base  currents  specified. 
The  distortion  in  the  distribution  of  current  across  the  array  decreases  as  the 
elements  approach  one-half  wavelength.  The  new  quasi  zeroth-order  theory 
shows  that  the  conventional  assumption  of  identical  current  distributions  is 
justified  only  for  the  case  of  half-wavelength  elements. 

The  new  array  solution  is  extended  to  arrays  with  elements  of  unequal 
length.  This  extension  follows  from  the  modification  of  representations  for 
the  coupling  terms  which  appear  in  the  integral  equation.  Two  separate  cases 
are  considered:  all  elements  driven  by  slice  generators,  and  arrays  with  at 
least  one  parasitic  element.  Numerical  calculations  of  the  two-element  para¬ 
sitic  couplet  confirm  the  well-known  phenomenon  that  the  main  radiation  beam 
may  be  directed  by  adjusting  the  length  of  the  parasitic  element. 

The  solution  for  the  array  with  elements  of  unequal  length  is  extended 
to  the  Yagi  array.  This  solution  reduces  to  a  simple  form.  Numerical  results 
are  presented  for  a  representative  three-element  Yagi  array. 
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6.  THE  PARTICULAR  SOLUTION  FOR  CURTAIN  ARRAYS 
WITH  HALF-WAVELENGTH  ELEMENTS 

The  general  functional  form  of  the  element  currents  given  in  2:52 
presents  some  difficulties  for  the  case  ph  =  — .  A  similar  difficulty  was 
encountered  by  King  [1,  p.  S449]  for  circular  arrays.  For  both  circular 
and  curtain  arrays  the  zeroth-order  expression  for  the  element  current 

TT 

assumes  the  indeterminate  form  0/0  when  ph  =  ~  .  This  indeterminate 
form  for  the  curtain  array  case  is  easily  demonstrated  by  examining  2:52 
given  by 


2tt 


r  F  (h)>F. 

CO  o  dr 


sinp(h  -  |2| ) 


+  j 


2Tr 


/<> 


F  (h)y. 
o  dr 


u  v^ 


(cospz  -  cosph) ' 
6:1 


From  the  form  of  the  S^functions  at  ph 


It  follows  that 


lim  =  lim  ' 

phvTT  ph>5 


10  0. 
0  10. 

C  .  .  . 


6,2 


The  indeterminate  form  of  the  element  current  follows  directly,  with  6:2  in 

6:1,  or 


,  2TTCospz 
po  dr 


.  2ttcos  pz 

^  0 
c'o  dr 


6:3 


The  circular  array  case  has  been  treated  by  King  by  differentiating 
with  respect  to  ph  both  the  numeiator  and  denominate,  of  the  one  dimensional 
form  of  6:1  .  .  This  method  is  too  cumberson  for  application  to  the  curtain 
array.  For  example,  this  technique  requires  that  each  element  of  the  matrix 
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be  differentiated  with  respect  to  ph.  A  simpler  method  is  presented 
for  the  special  case  ph  =  ^  beginning  with  the  original  form  of  the  integral 

Ct 

equations  given  by  2:1,  or 

4wi/  A  ,  (z)  =  Cm  .(z')K  (z,  z')dz' 

o  zk  J  /Ly  211  ki 

-h 

=  -j^(Cj^cospz  +  -Iv^j^sinpjzl)  6:4 

<po 

where 

e'jP^ki  ./ - 2 - T' 

=  — R—  ■ 

ki 

With  ph  =  -^  the  vector  potential  at  z  =  h  is  given  by  6:4,  or 

h  N 

“"‘'oV'i  =  j'  =  -Jr <“  "i^k'  • 

-h  i=l  ^  ° 


The  integral  equation  for  the  vector  potential  difference  W  (z)  =  A  ,{z)  -  A  ,{h) 

ZK  ZlC  ZxC 

is  obtained  by  subtracting  6:6  from  6:4.  This  vector  potential  difference 
vanishes  at  z  =  h  along  with  the  current.  The  rearranged  integr?.1  cquath  n  is 
given  by 


h  N 


47t  V  W  (z) 
o  zk 


I 

1=1 


I  .(z')K  ,  ,(z,  z')dz' 
,  zi  kid 


whe  re 


-j  cospz  -f  sinp|z|  -  1) 

ro  ok 


R 


ki 


^kih 


6:7 


6:8 


Note  that  the  form  of  the  right-hand  side  of  6:7  suggests  the  leading  terms 
for  the  current  given  by  King  [1,  p.  S449]  .  In  the  general  curtain  array  case 
the  integrals  in  the  rearranged  equation  Z:6  were  separated  into  two  groups 
depending  on  the  way  their  leading  terms  vary  as  a  function  of  z.  The  same 
separation  principle  is  applied  to  the  special  case  (3h  =  y  ,  with  one  group 

C4 

1- 

varying  approximately  as  does  M"  =  sin(3|zl-  1,  and  the  other  group  as  does 
F^;  =  cospz.  The  functional  forms  of  the  integrals  given  in  Z:12  -  Z:14  with 
Z:10  -  Z:ll  are  applicable  to  6:7.  The  current  in  each  element  is  separated 
into  two  parts  in  the  form 

I  .(z)  =  I  .(z)  +  I  .(z)  6:9 

Zl  Ul  VI 

where  by  definition  the  leading  term  of  the  two  pa  rtf  of  6:9  have  the  forms 

I  :  I  .(z)^  F  6:10 

VI  OZ  Ul  oz 


where 


oz 


=  sinp  |z(  - 1  ; 


F  =  cosBz 
oz 


When  6: 1C  is  substituted  in  6:7  the  integrals  that  occur  may  be  expressed 
as  follows  for  all  k  and  in  the  range  1  to  N. 


r  I  .(z)K,  .  ,(z,z')dz'  ,  I  .(z)  -  D,  (z) 

J  Ul  kid  kidu  Ul  kidu 

-h 


6:11 


II 


-h 
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I 

-h 


I  .(z)K,  .^_{z,z')dz'  .(z)  -  D,^_(z);  pb,  .  <  1 

VI  kidR  KidR  VI  kidR  ki 


6:13 


II 

i 


k 


As  in  the  general  case,  it  is  assumed  that  the  functions  are 

defined  such  that  the  differences  D  .(z)  are  small  enough  that  they  may  be 

lex 

neglected  in  a  solution  of  quasi-zeroth  order.  When  6:11  -  6:14  are  substi¬ 
tuted  in  6:7  with  only  the  leading  terms  retained  and  the  analysis  restricted 
to  the  case  where  pb  >1,  k  ^  i  (cf.  2:21  et  seq.  ) ,  then  the  following 
separation  into  the  two  groups  of  equations  may  be  carried  out: 


I  (z)  =  -j - c —  V  M 

vk'  '/  St'b  OK.  z 


u/n 

/o  dR 


6:15 


N 

I 

i=l 


I  ,(z)  =  -j^C,  F 
uk  V  k  z 


6:16 


whe  re 


=  4^^  6 

dR  kidR  ik 


6:17 


For  the  case  ph  -  ^  follows  directly  from  6:15  that  the  leading  term 
in  I^j^(z)  is  always  for  each  value  of  k.  Similarly,  it  follows  directly  from 
6:l6  that  the  leading  term  in  I  (z)  is  of  the  form  F  .  The  components  of  the 

Uiv  z 

element  current  must  then  be  of  the  form 


I  . 

VI 


(z) 


-j  A.M 


h 

y 

Z 


I  .(z) 

U1 


=  B.F 
1  z 


6:18 
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I  .(z)  =  -jA.M  +  B.F 
Zl  1  Z  1  z 


.  6:1< 


Since  is  real,  then  from  6:15  it  follows  that  A-  is  real  when  V  is  real 

dR  1  ok 

and  from  6:l6  that  B;  =  B._  +  jB_  is  complex, 

‘  iR  il 

The  form  of  the  constant  is  given  by  setting  z  =  0  in  6:4,  or 


h  N 


.  6:20 


The  constant  C  in  6:20  is  obtained  readily  by  substituting  6:19  in  6:20 
with  the  result 


where 


h 

(0)  =f  M*',.] 

kiv  J  z'l 


.K  (0,  z')dz' 

1  iCl 


^  .  (0) 

kiu 


n 


If  6:21,  6:18  and  6:19  are  substituted  in  6:15,  6:1^  the  result  in  matrix  form 
(cf.  2: 34  et  seq,  )  is 


- 1  [jA} 
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and 


where 


2it 


ok 


6:25 


kiv 


Sif^  (1 

kidv' 


6  )  +  i'F  6 

ik'  J  kidi  ik 


6:26 


and 


Also 


and 


Kidu 


(6) 

kidu  kiu 


[*] 

u 


* 

llu 

^  ■ 

1  2u  . 

INu 

21u 

22u 

*  ■  *2Nu 

* 

.  .  . 

Nlu 

N2u 

•  *  *NNu 

[i^]  = 

V 


llv  I2v 
21 V  22v 


.  h  h 

A 

.  Nlv  N2v 


iNv 

2Nv 

*NNv 


6:27 


6:28 


6:29 
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The  functions  are  given  by  6:11  -  6:14  and  6:22,  6:23,  or 


■■  -Re& (h.O)  -  E|^(h,0)j  -  ^^(h,  h)  -  E|^(h.  h)|| 

'^idl  =  '"'lb*”'"'  -  -  V”’”]} 

■^kidv  =  K'”-  ‘>]  ■  [v”’  >■'  ■  ^b'*-  ■'>1 


6:30 


6:31 


6:32 


'^kidu  =  P"  =  2 


6:33 


6:34 


^kiu<^)  =  C^(h.O), 


b  =  b 


ki 


6:35 


The  quasi-zeroth  order  form  for  the  element  current  with  (3h  =  —  is 
given  by  6:19  with  6:24,  6:25,  with  the  result 


/l  (z)l 

I  ^  1 


2tt 

■J  ^  vi/h 
<fo  dR 


|v^|(sinp)z|-  1)  -  j^^[#J-^[*J]|v^cos(3z  6:36 


The  coefficients  in  6:36  differ  from  those  given  by  King  [l,  p,S449]  for  the 
circular  array,  A  numerical  comparison  will  be  given  for  the  isolated 
element  withJ0u=  10.  The  values  of  the  SK  functions  for  this  case  are  given  by 
6:30  -  6:35,  and  it  follows  from  Mack  [2]  that 


6.9087 


6:37 
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=  0.  3842 

1  Idl 

6:38 

'4',  ,  ,  =  7.6637  -  j0.6331 

6:39 

=  -7.0754  +  jl.0932 

6:40 

.  (0)  =  8.  3517  -  jl.  8518 

1  lu 

6:41 

The  above  values  are  substituted  in  6:36  to  give  the  quasi  zeroth-order  current, 
or 


I  (z)  =  V 
z 


0.  OOZOcospz  -  j  I?.  0607cos(3z  +  2.  41 24(sinpfz|  -  l)j^xl0 


-3 

6:42 


where 


XL  =  10  ,  ph  =  I 


The  comparable  result  quoted  by  King  [1,  p.  S450']  obtained  from  the  general 
case  is  given  by 


[I^(z)] 


^  87cospz  -  j|6.  84cospz  +  2,  I7(sinp|z)  -  l)j^ 


X  10 


6:43 


The  input  admittance  from  6;42  is 


=  (10.0020  -  j4. 6480)x  lO'"^  mho 


.  6:44 


The  input  admittance  from  6:43  is 
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=  (sirs?  -  j4:'67)xl0‘^  mho 
The  second-order  theoretical  value  is 


Y 

o 


(9.  26  -  j4.  62)x  10 


mho 


6:45 


6:46 


Note  that  the  current  distribution  6:42  is  in  excellent  agreement  with  6:43 
as  shown  in  Fig.  6-1.  Also,  either  current  form  is  in  good  agreement  with 
the  approximate  second-order  curve  given  by  King  [3,  p.  Il6,  Fig.  22.9]. 

The  solutions  which  correspond  to  Cases  I  and  II  are  determined 
readily  from  6:24  and  6:25. 

Case  I. 

The  N  base  currents  are  specified,  or 

I  .(0)  =  I".(0)  +  jI'.(C)i,  1  =  1,  2,  3,  ...  N  .  6:47 

Zl  Zl  Zl 

The  N  base  voltages  are  required  which  maintain  these  base  currents. 

From  6:19  it  follows  that  at  z  =  0 

=  B.J,  t  j(B.,  +  A.)  ,  Ph=3  6:48 


or 


r'.(o) 

Zl 


6:49 


I'.(O)  = 

Zl  il 


+  A. 


6:50 
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With  6:25  in  6:24  together  with  6:49  and  6:50,  the  driving  voltages  are 
given  by 


Zir 


[4>J 


Ztt 


J  vj>h 
^o  dR 


-  -1 

[4  ]< 

I  (0)] 

u 

J 

6:51 


The  quasi  ze roth-orde r  element  current  is  given  by  6:19  with  the  A  coefficient 
given  by  6:25  and  the  B  coefficient  given  by  6:49  and  6:50  with  6:25. 

Case  II. 

The  N  base  potentials  k  =  1,2,  3,  .  .  .  N  are  specified,  and  the 

values  of  the  base  currents  I  (0)  are  easily  determined  from  6:36,  with  the 

z 

result 


Ztt 

J  ,•  \vh 
/o  dR 


Ztt 


fo  dR 


6:52 


The  analysis  of  the  three-element  endfire  array  with  quarter-wave 
spacing,  ph  =  -^  ,  jri,=  10,  and  specified  base  currents  follows  from  the  re¬ 
sults  for  Case  1  above.  The  analysis  parallels  that  given  in  Chapter  3  for  a 
similar  array  with  ph  =  ir.  It  is  expected  that  the  present  results  will  reduce 
approximately  to  the  conventional  case  where  the  distributions  are  identical. 
This  equivalence  is  inferred  from  the  similar  functional  form  of  and  and 
the  small  magnitude  of  the  ratio  (Aj^/B|)  for  the  isolated  element.  The  com¬ 
putations  for  the  driving-point  admittances  and  element  currents  are  given  in 
Appendix  VI,  and  summarized  below. 

The  driving-point  admittances  are  given  by  VI:  13,  or 


FIG.  6-1  QUASI  ZEROTH  ORDER  ELEMENT  CURRENT  FOR  )8h=T/2,fl  =  IO 
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^01 

=  (13.01 

-  j4.  v63)  X  10"^ 

mho 

^02 

=  (9. 393 

-  4.  939)  X  10'^ 

mho 

.  6:53 

^03 

=  (4. 322 

-  j3.  560)  X  lO"^ 

mho 

The  three-element  currents  are  given  with  respect  to  the  individual  driving 
voltages  in  VI:  14  and  with  respect  to  V^^  VI:  16.  They  are 


I  ,(z)  =  IC’^V 
zl 

1^2^^)  "  ■^V^2|9.  39  3cospz 

I  Az)  =  10‘\ 
z  5 


''oe 


{ 


Qj|^l  3.  Olcos  j3z  -  j  j?.  375cos(3z  +  2.  4l24{,sinpjz|  -  1^^ 

j  35 1  cos  pz  +  2.  41 24(sin p|z|  -  1  )j  ^  6: 
Q  2  322cos  pz  -  j  ^ .  972cos  pz  +  2.  41 24(sinp|z|  -  ^ 


54 


1 . 826(sin  pjz|-  1)  +  6,  7,65cos  p^  ”‘^1^’  ^  ~  1 )  -  9.  1 77cosp:^ 

9.  393cospz  ■  j  j^.  4124(sinp|z|  -  1)  +  7.  351cosp^^ 

1^2(2)  =  10  4759(sinp(zj  -  1)  -  9.  412cos  pz  -  ^  0.  9309(sinpjzl  -  1 )  +  10.  32cos  p:^ 

6:55^ 


The  element  currents  6:54  and  6:55  are  shown  in  Fig.  6-2.  The  element 

currents,  with  respect  to  the  individual  driving  voltages,  decrease  in  both  the 

real  and  imaginary  parts  from  element  one  to  element  three.  This  is  the  effect 

due  to  the  interaction  between  the  individual  element  currents  in  the  array. 

Note  that  the  general  shape  of  both  parts  of  the  element  current  are  similar 

in  all  three  elements  since  and  F  both  have  maxima  at  z  =  0  and  decrease 

z  z 

with  increasing  z.  When  the  element  currents  are  drawn  with  respect  to  V  , 

U  Cd 

the  distribution  of  the  current  on  all  three  elements  is  almost  identical.  Hence, 
the  quasi  zeroth-order  result  is  almost  identical  to  the  conventional  result 
for  this  special  case. 
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7.  EXAMPLES  OF  THE  GENERAL  THEORY: 


THE  INTERMEDIATE  CASE,  THE  THREE -ELEMENT  ENDFIRE  ARRAY 

WITH  Ph  =  AND  QUARTER- WAVELENGTH  SPACING 

The  three-element  array  with  ph  =  tt  has  been  covered  in  the  preced¬ 
ing  sections.  When  the  base  currents  are  specified  for  this  array,  the  shifted 
cosine  currents  are  automatically  determined.  However,  since  the  sinusoidal 
component  of  the  element  current  is  zero  at  the  driving  point,  it  is  not  affected 
by  the  specification  of  the  base  current  on  each  element.  The  situation  is 
somewhat  different  for  a  similar  array  with  ph  =  ,  since  both  the  sinusoidal 

and  shifted  cosinusoidal  components  of  the  element  contribute  to  the  base  cur¬ 
rent.  The  general  analysis  for  this  case  follows  a  similar  presentation  for 
the  full  wavelength  array. 

The  relationship  between  the  specified  base  currents  and  the  driving 
voltages  IS  given  by  2:dC.  The  numerical  calculations  are  carried  out  in 
Appendix  VII^  and  the  results  for  the  array  having  quarter-wavelength  spacing, 
.0.=  10,  and  the  main  beam  in  the  endfire  position,  are  summarized  below. 

The  driving-point  impedances  and  admittances  are 

=  357.0  +j293.7  ohms  ,  =  (1. 6706  -  jl. 3744  )x  10  ^  mho 

Z  =  374.  3  +  j402.  0  ohms  ,  Y  =  ( 1 .  3324  -  j  1 . 2402)xl  O’ ^  mho  7:1 

U  ^  V  C 

=  540.  8  +  j488.  5  ohms  ,  Y,  ,  =  { 1 . 01  83  -  jO.  9 1  V8)xl  C '  ^  mho 

0  3  L  J 

The  element  currents  are  given  with  respect  to  the  individual  driving  voltages 
from  VII:  1  3,  or 


10  -5  0  5  10  -10  -5  0  5  10 

MILLIAMPERES  PER  VOLT 


o 


FIGURE  6-2 
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I  ,{z)  =  V„,[0.9786F  -j(3.3219M  -O.SVr'SF  )]xlO 

zl  01  oz  oz  oz 


-3 


I  ,{z)  =  V^,[C.  7805F  -j(3.3219M  -  0.  6495F  )]xl0 

z2  02  oz  oz  oz 


-3 


7:2 


I  =  V^,[0.5965F  -j(3.3219M  -  0.  8371F  )lxl0‘ 

z3  03  oz  oz  oz 


where 


^oz  " 


F  =  cosRz  -  0.7071 
oz 


The  element  currents  are  given  with  respect  to  from  VII:  l6,  or 


7:3 


I  ,(z)  =  V„,[2.  7618M  -  b.  3669F  -  j(0.  3704M  -  0.  8798F  )lxl0 

zl  02  oz  oz  oz  oz 


-3 


=  \2l 


0.78C5F  -j(3.  3219M  -  0.6495F  )lxl0"^  7:4 

oz  oz  oz 


I  Jz)  =  ,[-4.  3902M  +  1.0378F  +j(0..3818M  -  0.  8844F  )]xl0 

z3  02  oz  oz  oz  oz 


-3 


The  above  results  are  easily  interpreted  by  investigating  the  behavior 

of  the  two  components  of  the  element  current,  F  and  M  .  From  7:2  note 

oz  oz 

that  the  real  part  of  decreases  and  the  imaginary  part  increases  from 

element  one  to  element  three.  This  result  is  identical  to  the  result  for  the 

same  array  with  Sh  =  it  (cf.  3:26  -  3:28).  As  before,  the  M  component  of 

oz 

the  element  current  is  proportional  to  the  driving  voltage.  Note  that  in 
7:2  the  imaginary  part  of  the  element  current  will  have  a  decreasing  amplitude 
at  z  =  0  from  element  one  to  element  three,  since  the  component  is  larg¬ 

er  than  the  F^^  component.  This  fact  explains  the  difference  between  the 
variation  of  the  admittances  for  the  case  ph  =  it  in  3:20  -  3:22  and  the 
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results  in  7:1.  That  is,  the  imaginary  part  of  the  admittance  decreases  for 
this  case  instead  of  increasing  as  in  the  case  (3h  =  it. 

The  element  currents  are  drawn  with  respect  to  the  individual  driving 
voltages  in  Fig.  7-la..  There  is  a  mild  decrease  in  the  real  and  imaginary 
parts  of  the  element  current  from  element  one  to  element  three.  To  compen¬ 
sate  for  this  effect  the  driving  voltages  must  increase  from  element  one  to 
element  three.  This  increase  will  affect  the  sinusoidal  current  component, 
since  it  is  directly  proportional  to  the  driving  voltage.  The  shifted  cosine 
component  of  the  current  will  also  be  affected.  However,  since  this  component 
•  ill  compared  to  the  s/TC..soi  lal  component,  it  will  not  be  as  effective 
in  changing  the  current  distributions  from  element  to  element.  The  change  in 
the  element  distributions  for  the  case  of  specified  base  currents  is  shown  in 
Fig.  7- lb,  where  the  currents  are  all  drawn  with  respect  to  Vq2*  Note  that 
the  variation  in  the  current  distributions  across  the  array  are  not  as  drastic 
as  those  given  for  the  same  array  with  ph  =  it  (cf.  Fig.  3-3). 


(b)  ELEMENT  CURRENTS  DRAWN  WITH  RESPECT  TO  Voz- 

FIG.  7-1  THREE  ELEMENT  ENDFIRE  ARRAY  BASE  CURRENTS  SPECIFIED 
fl  =  10,  SPACING) 
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8.  CURTAIN  ARRAYS  WITH  ELEMENTS  OF  UNEQUAL  LENGTH; 

THE  TWO- ELEMENT  COUPLET 
The  problem  of  determining  the  driving-point  impedance  for  two 
elements  of  arbitrary  length  and  spacing  has  been  carried  out  by  many 
authors  under  the  conventional  assumption  of  a  sinusoidal  distribution  on 
each  element.  A  history  of  this  problem  will  be  found  in  a  1951  article  by 
Chaney  [1].  Since  the  problem  of  two  skew  antennas  in  space  leads  to  many 
mathematical  and  experimental  difficulties,  this  discussion  will  be  limited 
to  parallel  elements  with  both  elements  driven.  The  solution  of  this  config¬ 
uration  leads  not  only  to  the  parasitic  couplet,  but  to  the  practically  impor¬ 
tant  case  of  the  Yagi  array.  The  couplet  with  elements  of  unequal  length 
is  shown  in  Fig.  8-1  for  elements  with  half  lengths  h^  and  h^  separated  by  a 
center  to  center  distance  of  b  .  With  the  aid  of  suitable  trigonometric  ap- 
proximations,  the  quasi  zeroth-order  solution  obtained  for  the  curtain  array 
may  be  extended  to  include  the  case  of  elements  of  unequal  length. 

The  two  integral  equations  for  the  currents  I  ,(z,)  and  I  -(z_)  in  the 

z  1  1  i.c  L 

elements  of  the  array  shown  in  Fig.  8-1  are  formally  written  down  from 
2:6  with  the  result 


-n  ^o  o  1 

8:1 


.  jt  1 'n ’"S  <“2^0.2  +  h^oz2' 


8:2 


“k  =  ^«ki  = 

i=l  i=l 


8:3 
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F  .  =  F  (z.)  -  F  (h.)  =  cosfiz.  -  cosSh.  8:4 

ozi  o  1  O  1  1  1 


M 


OZl 


G  {h.)F  (h.)  =  sinP{h.  -  (z.j) 
o  1  O  1  1  '  i' 


8:5 


For  the  array  with  elements  of  equal  length  the  currents  were  con-r 
veniently  separated  into  two  parts,  I  (z)  =  I  (z)  +  I  (z).  The  leading  term 

Z  V  tl 

of  the  I^(z)  distribution  was  found  to  be  proportional  to  =  sinp(h  -  jzj), 

and  the  leading  term  of  the  distribution  was  proportional  to 

F^^  =  cos(3z  -  cosph.  With  simple  sinusoidal  and  cosinusoidal  representa¬ 
tions  for  the  element  currents,  the  behavior  of  the  vector  potential  differences 
that  appear  on  the  left-hand  side  of  the  integral  equations  was  easily  calcula¬ 
ted.  The  vector  potential  differences  were  found  to  vary  approximately  in 

one  of  two  possible  ways  as  a  function  of  z.  One  group  varied  as  M  ,  and 

oz 

the  other  group  varied  as  forms  for  the  integrals  allowed  the 

original  integral  equation  to  be  separated  into  two  equations,  one  for  each 
distribution. 

The  formal  solution  for  the  unequal  element  case  parallels  the  general 
solution  of  the  equal  element  case  presented  in  Chapter  2.  However,  the  solu¬ 
tion  is.  complicated  by  the  behavior  of  certain  integrals  which  appear  in  8:1 
and  8:2.  These  integrals  behave  with  z  as  does  a  linear  combination  of  M 

oz 

and  The  behavior  of  all  the  integrals  on  the  left-hand  side  of  the  integral 

equations  is  first  investigated. 

The  behavior  of  the  integrals  in  8:1  and  8:2  with  Kernels  I  K  ,  ,  and 

zl  1  Id 

^z2^2c?  to  isolated  element,  was  investigated  in  Chapter  2. 

The  integrals  with  kernels  1  ,,  and  1  _K_,  ,  are  concerned  with  the  ele- 

zl  12d  z2  21d 

ments  of  unequal  length  and  must  be  investigated  for  all  values  of  pb.  A  good 
example  that  shows  the  major  effects  of  coupling  between  two  unequal  length 
elements  is  the  half-and-full-wavelength  couplet.  (N.  B.  Since  present  tables 


1  2 

FIG.  8-1  TWO  ELEMENT  ARRAY  WITH  ELEMENTS 
OF  UNEQUAL  LENGTH- 
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do  not  include  the  values  hi  =  ^/4  and  h  =  X/2  for  all  functions,  this  couplet 
is  approximated  in  some  cases  by  ^  =1.6  and  3.  2.  )  A  sinusoidal  cur¬ 

rent  is  assumed  on  each  element,  and  the  vector  potential  difference 
W  (z)  =  A  (z)  -  A  (h)  is  determined  for  both  gb  large  and  small.  Let 

Z  Z  Z  r  ^ 

h  =  X/4  and  h  =  X/2;  then  the  vector  potential  difference  along  element  one 

i  w 

due  to  a  current  sin(3(^  -  element  two  is  given  by 


'b'  2’ 


S.  z , 


Vi’  i> 


whe  re 


8:6 


0<z<-^,  b=b. 
1  4  K1 


8:7 


Both  the  real  and  imaginary  parts  of  8:6  are  shown  in  Fig.  8-2.  Note  that 
small  values  of  pb  are  not  included,  since  the  behavior  of  8:6  for  this  range 
is  discussed  in  Chapter  2.  Fig.  8-  2  shows  that  8:6  varies  approximately  as 
does  for  pb  >  1  in  both  the  real  and  imaginary  parts.  Hence,  a  simple 

trigonometric  approximation  serves  for  the  corresponding  integral  in  8:1. 

The  vector  potential  difference  along  element  two  due  to  a  current 
cospZj  on  element  one  is  given  by 


^ohi 


jcosp^K2^d<=‘2’^P'^^'l  =  V^l'^2)  -  VN'V 


8:8 


where 


0  ^  pz^  <  3.2  ,  phj  =  I 


8:9 
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Tables  of  the  function  C|^  are  not  presently  available  for  the  complete 

range  8:9>  and  the  calculation  of  8:8  is  based  on  a  representation  in  terms  of 

more  elementary  functions  [1,  p.  v7].  The  imaginary  part  of  8:8  is  shown  in 

Fig.  8-3  and  the  real  part  in  Fig.  8-4,  The  imaginary  part  of  8:8  behaves 

li.^e  F  ,  for  all  values  of  6b  1 .  However,  the  real  part  has  a  more  com- 
oz 

plicated  behavior  for  pb  >  3.  2.  The  reason  for  this  complicated  behavior  is 
found  in  the  curvature  of  the  wavefronts  near  element  two.  Since  at  the  sur¬ 
face  of  element  two  z^I  (z)  ^  E  ,  where  z^  is  the  internal  impedance  per  unit 

z  ^  z 

length,  some  insight  into  the  problem  is  gained  by  examining  E  ,(z  ■)  due  to  a 

Z  Cl  c 

current  cos^z  on  element  one.  The  electric  field  E  ^^^2^  been  derived  by 
King  [Z,  p.  528]  from  the  scalar  and  vector  potentials,  and  is  given  by 


^z2<^2) 


- 

J  4v 


-jsRih  ^-jpR2h 


R 


Ih 


R 


) 


2h 


8:10 


whe  re 


‘'ih  =  ^  '^2  ■  ‘’l’^ 
■*2h  = 


8:11 

8:12 

8:13 


The  magnitude  and  phase  of  8:10  are  shown  in  Fig.  8-5  with  respect 
to  the  reference  value  at  z^  =  0.  The  z  variation  of  this  magnitude  closely 
iesembles  that  of  a  shifted  cosine  for  all  values  of  pb  >1.2.  The  phase  of 
E  shows  a  variation  in  form  as  pb  is  decreased.  For  pb  >  u,  the  shape 
of  the  phase  front  is  ellipsoidal.  When  pb  is  decreased  to  pb  =  tt/2,  a 

i  o  1  C 

maximum  appears  at  pz^  =  it/ 2.  This  distortion  is  due  to  the  variation  in 


FIG.  8-2  vector  POTENTIAL  DIFFERENCE  ON  HALF  WAVELENGTH  ELEMENT  DUE 
TO  SINUSOIDAL  CURRENT  ON  FULL  WAVELENGTH  ELEMENT 
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FIG.  8-3  VECTOR  POTENTIAL  DIFFERENCE  ON  FULL  WAVELENGTH 
ELEMENT  DUE  TO  HALF  WAVELENGTH  ELEMENT  WITH 
CURRENT  COS.  fiZ  (IMAGINARY  PART) 
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FIG.  8-4  VECTOR  POTENTIAL  DIFFERENCE  ON  FULL  WAVELENGTH  ELEMENT 
DUE  TO  HALF  WAVELENGTH  ELEMENT  WITH  CURRENT  COS.  j9Z 
(REAL  PART) 
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FIG.  8-5  ELECTRIC  FIELD  ALONG  FULL  WAVELENGTH  ELEMENT  DUE  TO 
CURRENT  COS.  J3Z  ON  HALF  WAVELENGTH  CURRENT. 
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phase  velocity  along  the  ellipsoidal  wavefronts  produced  by  ^  lement  one. 

Q 

The  velocity  of  propagation  is  equal  to  v^  =  3x  10°  meters  per  second  along 
the  semi-major  axis  [Z,  p.  540],  However,  the  velocity  of  propagation  along 
the  semi-minor  axis  is  always  greater  than  or  equal  to  v^.  For  example,  with 
Ph  =  n/2  =  Pb  the  phase  velocity  of  the  wavefront  along  the  semi-major 

i.  1  u 

axis  is  1.  3  times  the  axial  phase  velocity.  The  distortion  of  the  surfaces  of 

Co.'  -t^nt  phase  prevents  the  approximation  of  the  element  current  1  _,(z_)  by  a 

z2  Z 

single  term  of  form  ^  However,  it  is  possible  to  approximate  8:8  by  a 

O  Zi  b 

suitable  combination  of  M  ,  and  F  Note  that  for  Bb,^  >  tt  a  one  term 

oz2  oz2  ^  12 

approximation  of  8:8  is  possible  (i.e.  in  terms  of  F  ). 

oz2 

A  simple  physical  explanation  serves  to  demonstrate  why  it  is  possible 
to  approximate  the  vector  potential  difference  8:8  by  a  combination  of  M 

oz2 

and  F  Element  two  is  illuminated  by  a  non-linear  wave  front  due  to  the 

oz2  ' 

distributed  currents  on  element  one.  Crudely,  the  response  of  element  t' m 
must  be  in  terms  of  its  fundamental  modes  which  have  sinusoidal  and  shifted 
cosinusoidal  forms.  When  h^  =  h^,  both  distributions  are  present,  but  the 
distortion  in  the  illuminating  field  over  the  length  of  the  receiving  element 
is  small.  Hence  little  of  the  received  power  is  placed  in  the  sinusoidal  mode. 
The  two  term  approximation  for  8:8  is  given  by 

W  (z  )  =  W*^sinp(h  -  |z  I )  +  W^(cos6z,  -  cosph_)  .  8:14 

u2  u  ^  Z  \  Z\  u  '  Z  2 

The  constants  and  W^are  determined  so  that  a  good  approximation  is 

given  for  W  (z  )  in  8:14.  Although  more  sophisticated  methods  are  avail- 

able  for  the  determination  of  W  and  W  (e.g,  a  least  mean  square  fit)  a 

u  u 

two  point  approximation  for  ^^(z^)  is  readily  accomplished.  Furthermore, 
the  two  point  approximation  includes  the  true  value  of  W^(0),  which  is  impor¬ 
tant  for  the  correct  value  of  input  power.  The  functions  M  _  and  F  ,  have 

oz2  oz2 

a  maximum  in  the  interval  0  <  z  <  h  .  Hence,  the  two  maximum  values 
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Tn  f 

of  W  (z^)  serve  to  determine  W  and  W  ,  or  for  Bh  >  tt/ 2  it  follows  that 
u  2  u  u  ■  ~ 

(N.B.  For  Bh  <  it/2,  M  _  has  a  maximum  at  z  =  0,  ) 

oz2 


W  (0)  = 
u  u 


sinph  +  W  (1  -  cosph  ) 


u 


8:15 


W^(h^  -  4)  = 


+  w|^(sinph^  -  cosph^) 


8:16 


Equations  8:15  and  8:l6  are  solved  simultaneously  for  and 

u  “ 


with  the  result 


^m  _  (sinph2  -  cos  ph?)W,i(0)  -  (1  -  cosph2)W„(h2  -  h 
u  sinph  (sinph  -  cosph  )  -  (1  -  cosSh,) 


=  -W..(C)  +  sinph,w.,(h^  -  h 

u  sinph^(sinph^  -  cosph^)  -  (1  -  cosph^) 


.  8:18 


The  excellent  agreement  between  8:6  and  8:14,  with  the  values  given 
in  8:17  and  8:18,  is  shown  in  Fig.  8-6.  The  approximation  8:14  is  quite 
general,  and  reduces  to  the  normal  one  term  approximation  for  8:6,  for 

<  1  and  >  rr.  Before  the  solution  for  the  original  integral  equa¬ 

tions  is  obtained,  the  behavior  of  the  integrals 


z'|)Kkid(z.  z')dz' 


8:  IV 


and 


(cospz  -  cosph)Kk^j(z,  z')dz' 


-h 


8:20 
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FIG.  8-6  COMPARISON  OF  TWO  TERM  APPROXIMATION  FOR  VECTOR 


POTENTIAL  SHOWN  IN  FIG.  8-5  WITH  TRUE  VALUE 
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is  given  in  Table  8-1  as  a  function  of  |3b.^"  An  inspection  of  Table  8-1  shows 
that  the  only  complication  arises  when  h^  >  h^^  for  pb  >1.  In  all  other  cases 
the  behavior  of  the  two  integrals  is  identical  to  the  equal  element  case. 


TABLE  8-1 


Integral 

Behavior 

Pb 

<  1 

Pb  > 

1 

Re 

Im 

^2 

<  hi 

h2  > 

hi 

Re 

Im 

Re 

Im 

M 

oz 

F 

OZ 

F 

oz 

F 

oz 

W^M 

V  oz. 

F 

oz 

F 

oz 

F 

oz 

F 

oz 

F 

oz 

W^M 
u  oz 

F 

oz 

-h 

+W^F 
u  oz 

The  solution  for  the  unequal  element  case  is  initiated  by  separating  the 
current  in  each  element  into  two  parts  of  the  form 

I  .(z)  =  I  .(z)  +  I  .(z)  8:21 

Zl  VI  U1 

where  by  definition  the  leading  terms  of  the  two  parts  of  8:21  have  the  forms 

I  .(z.)ro  M  .  ,  I  .(z.)  /\y  F  .  .  8:22 

VI  1  OZl  U1  1  OZl 

When  8:21  is  substituted  in  8:1  and  8:2,  the  following  integrals  occur 
(cf.  Table  8-1)  for  k  =  1,  2  and  i  =  1,  2  (N.  B,  The  higher  order  difference 
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tei’ms  have  been  neglected.) 


..h 


-h" 


r'  A 

\  .1  .(z;)K  .^^(z  ,z.')dz.'  -  'i'.  ,  (z,  ),  po,  .  <  1 

J  vi'  i'  KidR  K  1  1  kidR'A  vk  k' 

iC 


ki 


h 


-h 


[  I  .(z.')K  (z  .z.')dz;  =  ,  (z,  )(j^)  .  po,  .  <  1 

J  VI  1  Kidl  1  1  ivial  uk  k  Bj.  ^  ' 


ki 


i 


I  .(z.')K,  ,  ,(z,  .  z.')dz,'  ,  I  (z,  ),  h,  :^h. 

ai  1  kid  k  1  1  Kidu  k  k  i 


8:23 


8:24 


8:25 


i" 


I  .(z.')K  (z,  ,z.')dz.'  = 

VI  1  kid  k  1  1 


piD 


ki 


>  1 


8:26 


,h 


f  I  .(z.')K  .  _(z,';  z:)dz.'  =  4^""  (^)I  ,  +  4'^  .  {fi)I  .  .  h  >h.  8:27 

,J  ui  I  KidR  k  1  1  ivXdu  jA,  vk  iN.idu  E,  uk  i 

-h^  k  k 


t  ,  ,  h  >  h.  8:28 
uk  k  1 


B  * 

I  .(z.')K  dz  ,  z.')dz.'  ^  Inn4^  {•^)I  ,  (z  )  ,  h  >  h. 
ui  1  kidr  k  1  1  .vidu  Bj^.  uk  k  k  i 


8:29 


f  I  .(z.')K  (z  ,z.')dz.'  (i^)I  ,  ,  h,  >  h. 

,  J  VI  1  kidl  k  1  1  Kidv  B,  uk  .  k  i 

-h  k 


8:30 


where  K  .  and  K,  ,  are  defined  by  2:10,  or 
i.idR  kidl 


K,  .  j  =  K,  +  jK, 
kid  kidR  kidl 
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and  where  Im  and  Re  are  the  imaginary  and  real  parts,  respectively.  Also 

for  (3h  >  X./4 

ic 

^  (sinph|,  -  co3phk)Re^tt^^fiv(0)  -  (1  -  cos phk)Re ^idv(h^  - 
kidv  sinph  (sinSh  -  cos^h  )  -  (1  -  cosflh  ) 

As.  is.  rC 

"^idv^\  '  4^  _  8:32 

sin(3hj^(sinph  -  cosph^  )  -  (1  -  cosph^  ) 


The  constants  x.  , 
kidu 


and 


kidu 


are  given  by  8:31  and  8:32  with  ReSt^  .  , 

kidv 


replaced  by  Re'4f  .  ,  . 

kidu 

When  the  integrals  on  the  left-hand  side  of  8:1  and  8:2  are  replaced  by 
their  approximate  forms  from  8:23  -  8:3C,  and  it  is  assumed  that  the  elements 
are  spaced  greater  than  pb  =  1,  the  following  separation  into  two  groups  of 
equations  may  be  carried  out: 


=  J 


2Tr 


'fF 

o 


■V  .M 


TldR  vT  '  V  F  (h,)  or  ozl 


+  ( Aj)  ^±L.) 

**11  J..'  A  '  f 


22dR  2]dv'A-'  21du'jA_  '  v2  )  YF(h,)  02  oz2 
>—  2  2-1  3o  O  2 


8:33 


lldFB^^  12dv'  12du'B^' 


Iui(2)  =J 


4Tr 


■U,F 


•J^^F_(h,)  1  ozl 


o'“l' 


8:34 

“  V  F  (hj'^/kz 


At  this  point  the  simplicity  of  the  quasi  zeroth-order  solution  is  clearly 
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demonstrated,  since  it  follows  directly  from  8:33  that  the  leading  term  in 


the  I  ,  (z)  is  M  ,  for  each  value  of  k.  In  the  same  manner,  it  follows  directly 
vk  ozk 

from  8:34  that  the  leading  term  in  I  (z)  due  to  all  contributions  is  of  the 

“  Ulv 

form  F  .  It  is  now  possible  to  set 
ozk 


I  .(z)  =  jA.M 

VI  1  ozi 


B.F 
1  ozi 


8:35 


or 


I  .(z)  =  jA.M  .  +  B.F  •  8:36 

Zl  1  OZI  1  OZI 


Since  is  real,  it  follows  from  8:33  that  A,  is  real  when  V„,  is 

iidR  1  01 

real  and  A^  is  in  general  complex.  From  8:34  it  follows  that  in  general,  B^ 
and  are  complex. 

The  constant  U  is  obtained  readily  by  substituting  8:36  in  8:3,  with 
K 

the  result 


^k- 


-j/r. 

47T 


i 

i=l 


(h.  )  +  B.X-  (h,  ) 

V  k  1  kiu  k 


whe  re  k  =  1,2  and 


(h,  )  =  Cm  ,,  K.  .(h,  ,z')dz' 
Kiv  k'  oz'k  ki'  k 


-h 


8:37 


8:38 


8:39 


When  8:37,  8:35,  and  8:36  are  substituted  in  8:33  and  8:34,  the  result  is 
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Zn 


o'  r 

.m 


S  8:40 


•’^z'^ZdR  ^  •’^l^Zldv  ■'■  ®1^21du  F  (h  )  ^02 


X  ^  ( 

o 


[SlduW  -'^2u"h]®2=R.v<V  -i'Sldl"'o<>^j'j^l 

+  [Vh>  -  Wo<M«^ 

[l-^dv  "  j'""^du'^o<‘’2'  -  "|^u''o<V  ■'^u"'J^2  = 

[■^ivly  -  ‘'^dv+  i""'^2,dv'^o'>^  '  [*'2A-'y  - 


>8:41 


Note  that  the  physical  situation  is  somewhat  different  from  the  equal 
element  case  since  both  the  A  and  B  coefficients  contribute  to  the  sinusoidal 
distribution  of  current  on  element  two. 

Equation  8:41  may  be  written  in  matrix  form  with  the  following 
substitutions: 


’^iduVN>-'*'udu<N' 

VI/.,  F  (h  (h  ) 

12du  o  1  1 2u  1 

= 

f 

_!_*21du+i'"’'^21du»^o<’'2'  -  '21ul'’2' 

S/,,  F  (hJ-vi/_  (h  ) 
22du  o  2  22du  2_ 

'n2v<N'-^12duV>'ll' 

l«“l  = 

V 

>v<''2'-''^lldv*j''"'nidv> 

'n2v(''2'-J'n2dI^<''2|_ 
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With  8:42  and  8:43  in  8:41  it  follows  that 


Equation  8:44  resembles  the  result  2:39  for  the  equal  element  case. 

However,  the  relation  between  the  driving  voltages  and  the  A  coefficients  is 

more  complicated.  The  A  coefficients  are  related  to  the  driving  voltages  by 

introducing  8:40  into  8:44,  that  is,  by  suhratituting  the  value  of  A  from  8:40 

given  in  terms  of  B,  and  V  into  8:44.  The  value  of  A  is  given  by  8:40,  or 

1  o 


F  (h) 
o 


-It.l  IV  B 


where 


(♦.! 


'i'. 


lldR 

m 


12dv 


0 

22dR 


l*cl 


C 

k^m 

12du 


0 

0 


8:45 


8:46 


With  8:45  in  8:44,  the  B  coefficients  are  given  in  terms  of  the  driving  voltage 
by  the  matrix  equation 


r 


8:47 
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or 


.2Tr 


l»“l 

u 


“1-1  -if  "1 

]  "  i-  ^-Q  i 

.j  ^F^(h)r 


8:48 


The  solution  of  the  array  problems  which  correspond  to  Cases  I  and  II  for 
the  array  of  unequal  lengths  may  be  obtained  in  the  same  manner  as  shown 
in  Chapter  2. 
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9.  CURTAIN  ARRAYS  WITH  ELEMENTS  OF  UNEQUAL  LENGTH; 
THE  TWO-ELEMENT  (COUPLET,  SPECIAL  CASES 
The  special  form  for  the  element  currents  in  a  curtain  array  for  the 
special  case  ph  =  ir/Z  is  given  in  Chapter  6.  Thus,  from  6:9  and  6:10 


I  .  =  -jA  M^  +  B.F 

Z1  i  OZ  1  oz 


9:1 


where 


oz 


sin  Pjzj  -  1  , 


F  =  cosQz 
oz 


9:2 


The  special  form  of  the  element  current  in  9:1  was  obtained  by  examin¬ 
ing  the  original  integral  equation  for  the  currents,  2:1,  with  the  value  of  ph 
set  equal  to  rr/Z.  Since  curtain  arrays  with  at  least  one  element  of  half  length 
ph  =  Tr/2  are  quite  common  in  practise,  it  is  important  to  investigate  this 
case  in  detail.  Consider  the  two-element  array  shown  in  Fig.  9-1,  with  one 
element  of  half  length  ph  ^  -nil  and  another  with  ph  =  tt/Z  driven  by  potential 

*  b 

differences  and  The  two  integral  equations  for  this  couplet  may  be 

written  directly  from  2:6  and  6:7.  Thus 


“l  -^2  o  V 


9:3 


^1  -ho  02 
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where 


U, 


)  -jfa  \i.(z.')K,.(h,,2.')dz,' 
4t:  J  zi.  1  li  1  1  1 


9:5 


)K,.(0,2.')dz.' 
Zi  1  1 


9:6 


g'jP^ki  e 


9:7 


ki 


kih 


Before  the  behavior  of  the  various  integrals  in  /:3  and  9:4  is  inves¬ 
tigated.  it  is  useful  to  give  a  qualitative  description  of  the  self  and  mutual 
terms.  With  no  couplingi  the  current  on  element  one  is  given  by 

I  ,(z)  =  jA.M  ,  +  B  F  and  I  (z)  =  -jA  _  +  B  F  for  element 
zl  1  ozl  1  ozl  z2  2  oz2  2  oz2 

two.  When  these  two  driven  elements  are  placed  in  proximity,  currents 
are  induced  on  each  element  from  the  opposite  member.  Two  cases  may  be 
distinguished:  (3h^  <  Tr/2  and  ph^  >  v/Z.  For  ph^  <  tt/2  a  current  propor¬ 
tional  to  F  ,  is  produced  due  to  both  parts  of  the  current  on  element  two. 
ozl 

A  more  complicated  behavior  is  exhibited  for  the  current  induced  on  element 

two  due  to  the  current  on  element  one.  Due  to  the  distorted  wavefront  which 

illuminates  element  two,  an  induced  current  is  produced  which  contains  a  term 

proportional  to  F  plus  a  term  proportional  to  M  .  For  the  case  ph  >  it/2 

the  behavior  of  the  individual  elements  is  reversed.  The  induced  current  on 

element  one  contains  terms  which  are  proportional  to  F  ,  and  M  How- 

ozl  ozl 

ever,  on  element  two,  a  current  proportional  to  F  is  produced  by  I  (z). 

OZ  ^  Z  1 

For  both  cases,  it  is  the  larger  element  which  exhibits  a  complex  behavior 
while  the  smaller  acts  like  a  receiving  element  in  a  uniform  field. 
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The  current  on  each  element  is  first  separated  into  two  parts  of  the 


form 


I  .(z)  =  I  ,(z);+  I  .(z)  i  =  1,  2 

Zl  VI  Ul 


9:8 


where  by  definition  the  leading  terms  of  the  two  parts  have  the  forms 


I  (z)  M 
vi  ozl 


9:9 


‘u2<^'  ^oz2 


.  9:10 


When  the  above  current  forms  are  substituted  in  9:3  and  9:4,  the  fol¬ 
lowing  integrals  occur  and  may  be  expressed  as  shown  in  Chapter  8 
(c.f.  8:23  -  8:32),  where  the  higher  order  differences  have  been  neglected: 


(  ^vl^^l^^lldR^^l’ ^P^^l  ^  ^lldR^vl^^l^  ’  •^*^11^  ^ 


9:11 


rr  A 


9:12 


r 


J  'v2<*2>‘S2d<'l’“2W 


'^l''2dv<'^‘ulln''  ”1  "  *’2 


Re 


9:13 


9:14 


hj  >  h2 


V. 
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r 


^12du’'B^*'ul*“l*  ’  ’’l  ^  '’2 


9:15 


hi  >  h^ 


ill<^''K21clV^''‘*^l'  = 


Re  I 


/iAj 


9:16 


'^12dv^B2  ’  ^1  ^  ^2 


Re 


'*'l2du^B^^^u2^^2^  *  ^1  ^  ^2 


,B 


9:17 


->  ..h 


-h‘ 


J^v2^^2^^22dR^^2’ ^2^^^2  ""  '^22dR^v2^^2^ '  *^*^22  ^ 


9:18 


-h' 


j'v2<^z'‘'22dl'=‘2'  ^^'“^2  -  ’^22dl‘^'‘u2'^2'  ’  ^^22  "  * 


v:19 


J 


^u2^^2^^22d^^2’  ’^''i}^'^'l  '^22du^u2^^2^ 


9:20 


where 


_  0.7071Re4/i?.H,.(0)  +  0.  2929Re ,Hn(t 

12du  -0.4142 
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sf4  =  -Re^i2du(0)  -  Re'»i2du(f)  o.^z 

12du  -0.4142  •  ^ 

For'4','?,  and'4',^,,  ,  replace  Sk  by  '4'.  _  ,  in  9:21  and  9:22.  The  two  term 

12dv  12dv  ^  12du  12dv 

approximation  for  the  real  part  of  '4^2du  obtained  from  a  two  point  approxi¬ 
mation  at  =  0  and  pz2  =  n:/4  (c.f.  8:14  et  seq.  ). 

The  integrals  on  the  left-hand  side  of  9:3  and  9:4  are  replaced  by  the 
approximate  forms  given  in  9:11  -  9:20,  and  for  pb^^  ^  ^  following  separa¬ 
tion  into  two  groups  of  equations  may  be  carried  out  for  h^  <  h2.  (N.  B.  Higher 

order  terms  have  been  neglected.  ) 


2ti 


'^IdR^vl^^^  V  F  (h 
po  o  1 


^^l^ozl 


9:23 


22dR 


21dv'  A  2' 


+  SK 


m 


B 


21du'-jA2'l  v2 


}v2 


(z) 


h 

oz2 


9:24 


Idl  '  B  '  '  ’ 


i  •:  J  po  o  1 

[^'^Idv'^-’^'^'^ldv^  ^Idu  ^  B2^'^2dl  '^'^2dv^li2^^^  ""  ~^^VozZ 


9:25 


9:26 


It  follows  directly  from  9:23  that  the  leading  term  for  I^^(z)  is  pro¬ 
portional  to  3,nd  from  9:24  that  the  leading  term  for  I^^(z)  is  propor¬ 

tional  to  ^^^2'  ^i'^i^l^.rly ,  it  follows  from  9:25  that  the  leading  term  for 

I  ,(z)  is  proportional  to  F  and  from  9:26  that  I  ^(z)  is  proportional  to 
ul  ozl  u2 

F  Hence,  it  is  now  possible  to  set 

oz2 
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vl<“>  =  j^l“ozr 

I  (z)  =  B  F 
uT  '  1  ozl 

9:27 

V2<*>  =  ■j^z'^^zZ  ■ 

I  -(z)  =  B  F  , 
u2  2  oz2 

9:28 

=  -j^Z^o^zZ  *  ^Z^'ozZ 


9:29 

9:30 


Since  is  real,  it  follows  from  9:23  that  A,  is  real  when  V„,  is  real. 

lldR  1  01 

However,  9:24  shows  that  is  in  general  complex.  From  9:25  and  9:26, 

it  follows  that  in  general,  and  are  complex. 

The  constants  and  are  given  by  substituting  the  element  currents 
given  by  9:29  and  9:30  into  9:5  and  9:6,  with  the  result 


U, 


-AFj 

Mit  ^ 


JAi'i'nvO',) 


- 

•’4Tr 


V;31 


;  .A 

2  •^4-11 


jA^4;^lv(0)4B^HK^lu(0) 


+  j^ 


-jA2^,2v<°)"®2'^22u(°) 


9:32 


With  9:31  and  9:32  substituted  in  9:23  -  9:26,  the  result  is 


i  A  ^  -  i  _ V 

•'  1  lldR  V  F  (h,)  01 

Ao  o  1 

•’^I'^Zldv  ®l'*'^du  ■  ■^^z'^tzdR 


-jI'v 


Jc 


02 


9:33 
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po<VSlldu-'Slu«H|®l  +  [o<V'^u-'^12u'vj®2  =  [W-^''o<V'^ld^'j^> 

9:2 

Equations  9:33  and  9:34  are  easily  written  in  matrix  form  as  follows: 

|t,l|iA|  =  j^l^}  9:3 

and 


where 


IV  = 


‘P  0 

lldR 


q/m  q/ 

-  21dv  22dR 


[♦cl  - 

^21dv 


0 


9:38 


SR12 


-83- 


[^] 

u-* 


1<1  = 


F  -'F,  (h,)' 

o  1  12du  12u  1 

K2du*J'"'^21du>-'*'21u<“' 

■^22du-'^Z2u<"' 

Klv<N'-j"'o"'.>^lldI 

•^!'2v'>‘l'-'’o'S'-*!'2dv] 

9:39 


-  j-*-22dI 


9:40 


The  result  comparable  to  9:35  and  9:36  for  the  case  h  h  has  the 

1  w 

following  values  for  the  component  matrices 


1*1 


1  IdR 


p,m 

‘12dv 


4^ 

22dR 


9:41 


[^,1 


^tn 

12du 


9:42 


I*,,  I 


F  (h,)-d^,  ,  -  -F;,  (h  ) 

o  1  1  Idu  1  lu  1 


'i'  -  4-  (0) 

12du  21u'  ' 


'F  -  (0) 

22du  22u'  ' 


9:43 


q-44 
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The  matrix  equations  for  the  couplet  given  by  9:35  and  9:36  are 
coupled.  That  is,  both  the  A  and  B  current  coefficients  appear  in  both 
equations.  When  both  elements  are  equal  in  length,  the  coupling  disappears 
and  the  sinusoidal  current  coefficient  jA^  again  is  proportional  to  the  driving 
voltage.  For  example,  one  matrix  equation  is  given  for  the  B  coefficients 
in  terms  of  the  driving  voltage  by  combining  9:35  and  9:36,  or 


.  2tt' 

%■ 


9:45 


Consider  the  practical  case  where  it  is  desired  to  constrain  the  driving 
voltages  and  determine  the  corresponding  element  currents.  The  B  coeffi¬ 
cients  of  the  element  cur  rents  are  given  by  9:45  in  terms  of  the  two  driving 
voltages  and  sinusoidal  current  coefficients  are  given  by  9:35 

with  9:45.  The  current  on  each  element  follows  directly  from  9:1. 
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10.  CURTAIN  ARRAYS  WITH  ELEMENTS  OF  UNEQUAL  LENGTH: 

THE  PARASITIC  CASE,  EXAMPLES  OF  THE  THEORY 

The  preceding  analysis  has  been  concerned  with  all  elements  driven 
by  separate  generators.  When  the  elements  are  equal  in  length,  the  sinusoidal 
current  on  each  element  is  proportional  to  the  driving  voltage.  This  is  the 
result  of  the  quasi  zeroth-order  approximation.  The  situation  is  somewhat 
changed  for  the  case  of  unequal  length  elements.  The  more  complicated  be¬ 
havior  of  the  coupling  terms  in  the  integral  equations  requires  sinusoidal  cur¬ 
rent  as  well  as  shifted-cosinusoidal  currents  for  a  good  approximation  to  the 
actual  behavior  of  the  integrals.  The  fundamental  approach  used  in  the  quasi 
zeroth-order  solution  is  a  separate  iteration  for  both  the  sinusoidal  and  cosinu¬ 
soidal  currents  in  the  array.  Hence,  for  the  array  with  elements  of  unequal 
length,  the  sinusoidal  currents  in  the  coupling  terms  are  determined  by 
suitably  modifying  the  equations  which  normally  account  for  self-induced 
sinusoidal  currenta  The  modification  involves  an  adjustment  of  the  driving 
voltages,  which  in  turn  compensates  for  the  effect  of  the  sinusoidal  currents 
in  the  coupling  terms. 

The  quasi  zeroth-order  theory  implies  that  the  sinusoidal  current  on 

an  element  in  an  array  vanishes  if  the  driving  voltage  is  zero.  This  result 

is  easily  seen  from  8:33  or  9:23  -  9:24  with  =  0,  i  =  1,  2.  This  is  only  an 

approximation,  since  the  simplification  is  made  that  certain  integrals  in  2:1 

vary  only  with  This  leads  to  a  discontinuity  in  the  approximate  form  of 

the  integral  at  z  =  0 .  whereas  this  discontinuity  does  not  exist  in  the  original 

integral.  However,  since  the  trigonometric  approximations  for  the  integrals 

are  made  at  the  point  of  maximum  current,  the  approximation  should  be  quite 

good  except  near  z  =  0.  One  method  of  checking  the  validity  of  assuming  M  =  0 

oek 

when  is  to  compare  the  approximate  result  with  a  more  rigorous  solution 

given  by  King  [1]  using  a  symmetrical  component  solution  with  N  =  2.  The 
solution  of  the  couplet  array  with  one  parasitic  element  is  given  as  a  super- 
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position  of  two  separate  couplet  problems,  each  with  prescribed  driving 
voltages.  First,  both  elements  are  driven  in  phase  with  voltages 
V  =  V  /2  =  V  ,  and  then  driven  in  phase  reversal  with  ,  =  V-/2  =  -V  . 

The  superposition  of  the  currents  for  these  two  cases  yields  the  required 
couplet  with  and  “  0.  An  interesting  case  to  compare  the 

approximate  and  rigorous  solutions  for  the  two-element  couplet  is  6h  =  3it/4 
and  Pb  =  ■rr/4.  The  calculations  are  given  in  Appendix  VIII.  The  result  which 
follows  from  the  quasi  zeroth-order  theory,  where  it  is  assumed  that 
zero,  is  given  in  VIII;i5,  or 

=  lo'\,  10.8030F^^-j(3.3219M_^^-  0,7282F^^)l 
-3 

I  ^(z)  =  10  V^,  [0.4732F  +  jO.  l604F  ]  , 

z2  01  oz  ■'  oz^ 


The  corresponding  result  based  on  a  rigorous  solution  in  terms  of  symmet¬ 
rical  components  is  given  by  VIII: 21,  or 


I  ,(z)  =  10'\„,  [0.9368F  -j{3.  3219M  -  0.6789F  )] 

zl  01  oz  oz  oz 

I  ^(z)  =  lO-^V^,  [0.6067F  +j0.1075F  1 

z2  01  oz  oz 


10:2 


The  element  currents  10:1  and  10:2  are  illustrated  in  Fig.  10-1  for  both 
the  real  and  imaginary  parts. 

The  error  at  ph  =  3'tT/4  should  be  nearly  maximum.  This  may  be 
deduced  by  examining  the  cases  6h  =  rr/Z  and  ph  =  tt.  At  Bh  =  tt/Z  the  element 


current  is  almost  entirely  composed  of  one  component,  the  F  distribution. 

oz 

At  ph  =  TT  the  parasitic  coupling  is  extremely  weak,  as  shown  by|King  1,  .p  S45'^ 
Hence,  for  both  these  extreme  cases  the  error  should  be  less  than  for  the 


case  ph  =  3Tr/4. 


L. 


^b=T 


MILLIAMPERES  PER  VOLT 


FIG.  lO-l  TWO  ELEMENT  ARRAY  WITH  /9h=3ir/4,  j8b=ir/4, 
ONE  ELEMENT  PARISITIC,  ELEMENT  CURRENTS. 


j8h=3ir/4 
V, 


MO 


)8h=3Tr/4 

V20=0 


■^10  O 

T 


0^J0  + 
z 


v,oi 

T 


0 


SR12 


-87- 


A  well-known  phenomenon  discussed  in  the  literature  [2]  is  that  for  a 
fixed  spacing  it  is  possible  to  direct  the  main  beam  in  either  direction  with 
respect  to  the  parasite.  This  directional  property  is  accomplished  by  adjus¬ 
ting  the  length  of  the  parasitic  element,  with  the  driven  element  fixed  near 
ph  =  Tv/Z.  This  effect  has  been  verified  experimentally  as  well  as  theoretically 
by  an  approximate  solution  due  to  King  [l].  However,  the  approximate  solu¬ 
tion  is  restricted  to  small  changes  in  the  length  of  the  parasitic  element  near 
ph  =  tt/2.  The  present  theory  is  not  restricted  with  respect  to  the  length  of 
the  parasitic  element. 

The  first  example  is  the  couplet  with  ph^  =  1.4,  6h^  =  tt/2,  =  0, 

V  .  The  complete  calculations  for  the  element  currents  are  carried  out  in 

U  Ct 

Appendix  IX,  with  the  result  given  by  IX: 28  and  IX: 29,  or 

I,  (z)  =  l0"^V„.,(-20.  78  +  17..  81)F  10:3 

zl  02  ozl 

I  .,(z)  =  lO'^V,  ,[(-0.6554  -.32.468)1^1^  -+  (14.  06  -  j25.  39)F  .  10:4 

ZiM  0  OZw  OZm 

The  element  currents  are  shown  in  Fig.  10-2,  drawn  with  respect  to 

The  driving-point  admittance  and  impedance  of  element  two  is  given  by  10:4, 

evaluated  at  z  =0. 

=  Y  =  10'^(14.72  -  j22.92)  mho  10:5 

^^02 

Z  =  TP —  =  {19.  84  +  j30.  89)  ohms  .  10:6 

^  ^02 

The  radiation  field  in  the  equatorial  plane  is  given  by 
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eJ(^) 


K;lC, 


10:7 


where 


C,  =  B  G  {-T,  1.4)  =  0.7475B, 
1  1  m  Z  1 


C,  =  -jA  H  +  b  G  =  -jA  (-0.  5708)  +  B, 

Z  I  xlx  Z  Z  1  ^ 

.  TT,  ainQ  -  sin(-^co36) 
m'  ’2^  ~  sinGcosO 


>  10:8 


With  10:3,  10:4  and  10:8  in  10:7  the  numerical  value  of  the  radiation  field 

reduces  to 


Eg(4>)  =  K[(-l.  553  +  jl.  331)  +  (1. 443  -  j2.  398)e^^°®‘*  ]  10:9 


where 


(3b  =  1 


Q 


10:10 


The  radiation  pattern  is  plotted  in  Fig.  10-3,  and  the  back  to  front  ratio  is 
given  by  ) 0:9»  or 


Ie(0)I 

|E(t0| 


ll. 245  +  jl. 2491 
1-2.791  -  ji.  179! 


0.  5822  'V  -4.  7  db 


.  10:11 


Note  that  for  this  case,  the  mam  beam  is  directed  toward  the  parasite.  The 
numerical  value  for  the  back  to  front  ratio  is  close  to  the  experimental  value 
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of  -6  db»  given  by  Starkey  and  Fitch  [3] . 

The  second  example  is  the  couplet  with  ph  =  1.8,  ph  =  tt/Z, 

X  Ld 

V  =  0,  and  The  calculations  for  the  element  currents  are  carried 

U  X  U  Ct 

out  in  Appendix  IX.  The  results  are 

I  ,(z)  .=  10’\„^(0.5701  +  j4.6794)F  ,  10:12 

zl  02  ozl 

I  Jz)  =  10"  V^-l-j2.4l24M^  ,  +  (7.178  -  j8. 865)F  ,]  ,  10:13 

The  element  currents  are  shown  in  Fig.  10-4,  drawn  with  respect  to 
Vq^i*  The  driving^oint  admittance  and  impedance  of  element  two  are  given 
by  10:13,  or 

=  10'^(7.  178  -  j6.453)  mho  10:14 

''02 

Z  =  77.05  +  j69.26  ohms  .  10:15 

u  ^ 

The  radiation  field  in  the  equatorial  plane  is  given  by 

E^{*)  =  K{C^  +  10:16 


whe  re 


C,  =  B  G  (^,1.8)  =  1.3837B,  10:17 

1  1  m  2  1 

''z  =  =  -iA;(-0.5708)  +  B^ 


^  10:18 
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With  10:12,  10:13,  10:17  and  10:18  in  10:l6,  the  numerical  value  of  the 
radiation  field  in  the  equatorial  plane  is  given  by 


Eg(*)  =  K[(0.  7888  +  j6.475)  +  (7.  178  -  j7.489)e-’‘^‘^°®*l  •  10:19 


The  front  to  back  ratio  is  easily  calculated  from  10:19: 


E(»  =  0)1  I1O.97  +  j8.469  I  _  14.21  _ 

|E(«=0)|  "  |-1.634  -  j3.6llj  “  3.965  "  ‘ 


.  10:20 


The  ratio  |E(0)|/  |E(Tr)|  in  10:20  corresponds  to  a  front  to  back  ratio 
of  11.  1  db.  Note  that  in  this  case  (Fig.  10-5),  the  main  beam  is  directed 
away  from  the  parasite.  The  front  to  back  ratio  of  11.  1  db  compares  favor¬ 
ably  with  a  value  of  11  db  given  by  McPetrie  and  Saxton  [4]  . 


-10  -5  0  5  10  -10  -5  0  5 

MILLIAMPERES  PER  VOLT  MILLIAMPERES  PER  VOLT 


FIG.  10-4  ELEMENT  CURRENTS  FOR  PARASITIC  COUPLET, 
i8h,*l.8,  1.4,  ;9b=  I 


SR12 


-gi¬ 


ll.  THE  YAGI  ARRAY 

The  general  Yagi  array  problem  is  considered  by  King  [l,  Ch.  VI, 

Sect.  5].  The  N  element  Yagi  array  normally  consists  of  one  driven  element, 
usually  a  half-wave  dipole,  one  reflector  of  half  lengh  Bh  >  tt/Z,  and  N-2 
directors  of  half  length  ph  <  ir/Z.  Since  the  arrsy  has  only  one  driven  element, 
it  has  an  advantage  over  the  conventional  endfire  array,  wherein  all  elements 
must  be  driven  in  the  correct  time  phase.  A  six-element  Yagi  array  is  shown 
in  Fig,  11-1. 

For  convenience  the  driven  element  is  chosen  to  be  of  half  length 
ph  =  ir/2,  and  the  reflector  {element  one)  of  half  length  ph  >  tt/2.  The  inte¬ 
gral  equations  for  this  array  easily  follow  from  the  results  for  the  couplets  of 
unequal  length 


4n 


i=l  'iu 


z,.z  )dz  =  V  U,F  , 

lid  1  i'  i  •’jf  F  (h, )  1  ozl 

C*o  o  1 ' 


^  (’l^.(z!)K,.  ,(z,,z.')dz!  =  -j|^V.-(-^F  ,  +  ,) 

i^  -hi  eo  02 

(So  o  3 


11:1 


11:2 


11:3 


i=l  -hj^ 


► 


li, 


hi 


Z1 


4tt 


U..F 


^oFoih^)  N  ozN 


11:4 


The  reflector  and  directors  in  the  Yagi  array  have  half-lengths  which 
are  near  ph  =  tt/Z  in  order  to  achieve  the  correct  directive  properties.  With 
this  in  mind,  it  is  possible  to  make  certain  simplifying  approximations  for 
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the  case  where  the  elements  are  spaced  a  distance  greater  than  pb  =  1.  For 
elements  of  half  lengths  near  ph  =  tt/Z,  it  is  assumed  that  all  induced  'currents 
are  proportional  to  That  is,  it  is  assumed  that  only  the  individual  ele¬ 
ment  produces  a  current  which  is  proportional  to  The  validity  of  this 

assumption  is  easily  seen  by  noting  that  near  ph  =  tt/2  the  and  distri¬ 

butions  are  similar  in  form.  Furthermore,  calculations  of  the  element 
currents  for  the  couplets  of  Chapter  10  show  that  the  corrections  for  not  assum¬ 
ing  certain  integrals  proportional  to  F^^  are  quite  small,  actually  of  the  order 
of  a  few  per  cent. 

An  example  of  the  order  of  magnitude  of  the  correction  term  is  given 

in  IX:27  and  IX;29  for  the  couplet  with  ph  =  1.4,  pli  =  tt/Z,  pb  =  1. 

The  reduction  of  the  original  integral  equations  for  the  Yagi  array 

11:1  -  11:4  to  a  solution  of  quasi  zeroth-order  follows  directly  from  the 

preceding  analysis  for  the  two-element  couplet  in  Chapters  8,  9  and  10.  The 

final  result  is  similar  to  9:35  and  9:36,  and  with  the  simplifying  assumption 

of  only  F  type  currents  for  the  coupling  terms,  the  element  currents  are 
oz 

given  by  the  solution  of  the  following  equations: 


11:6 


FIG.  Il-I  A  SIX  ELEMENT  YAGI  ARRAY 
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where 


llu 

12u 

.  . 

iNu 

f  ^ 

21u 

22u  ‘ 

•  *  2Nu 

i. 

Nlu 

N2u  ‘ 

‘  ■  NNu 

■^o'V-Sidu-W 

■  •  ^o^V'^u''^uV 

'^22du  ■'*^22u^^^  ■  ■  ■ 

■  ■  '^2Ndu"^2Nu^‘^^ 

■  ■  "'o<N>^Ndu-^Na"’3' 

/o‘V%ldu-1ilu<V 

■  ■ 

*  *  'W  V  ■ 

'^21v<^)  -^21dv 

-^zavl'^' -  •  ■  ■ 

’  • '^2Nv^^^ ‘"^Ndv 

-  "'o'N^ldv 

■%.v<W  - 

'WV  ■  ^o%^'^2iv  • 

■  •  "^W  V  ■ 

Note  that  the  formal  result  11:5  and  11:6  for  the  Yagi  array  is  iden¬ 
tical  to  that  for  the  N_element  curtain  array.  Also  note  that  since  only 
is  non-zero,  the  second  column  of  11:8  enters  into  the  calculations  for  the 
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element  currents. 

A  representative  three-element  Yagi  array  is  shown  in  Fig.  11-2, 
with  phj  =  1.8,  ph^  =  u/2,  ph^  =  1.4  and  pb  =  1  spacing.  The  calculations 
for  the  driving  impedance  and  element  currents  are  given  in  Appendix  XI. 

The  element  currents  are  given  by  XI:  9,  or 

I  ,{z)  =  10‘V^  f  (0.7355  +  j6.4421)F  ,]  11:9 

zl  02  ozl 

I  .(z)  =  10"^V^_[-j2.4124M^  _  +  (13.  19  -  j25.  55)F  1  11:10 

2  m  0  m  02  M  O  2  M 

I  ,(z)  =  10‘\^,  [  (-19.  36  +  jl4.  08)F  ]  11:11 

Z3  04  ozS 


where 


^o^z2  =  ^ 

11:12 

F  =  cosPz  -  cosPh, 

OZK  k 


The  driving-point  admittance  and  impedance  are  given  by  11:10 
evaluated  at  z  =  0.  Thus 

=  (13.  19  -  j23.  14)  X  lO"^  mho  11:13 

Z,_  =  18.59  +  j  32. 62  ohm  .  11:14 

02 


The  radiation  pattern  for  this  array  is  shown  in  Fig.  11-3.  The  front 
to  back  ratio  of  10.  1  db  is  less  than  optimum  for  this  array  and  is  nearly 
equal  to  the  corresponding  ratio  for  the  two-element  parasitic  couplet  with 


)3h,=  l.8  )0h2=Tr/2  iShjM.S 

a  =  10  FOR  ALL  ELEMENTS 


FIG.  11-2  THREE  ELEMENT  YAGI  ARRAY 


IK 
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Ph  =1.8,  ph  =  it/2.  Hence,  for  maximizing  the  front  to  back  ratio  of  the 

1  Ct 

three-element  Yagi,,  the  relative  positions  of  the  director  and  reflector  must 
be  changed. 

A  further  theoretical  simplification  is  available  for  the  Yagi  array, 

since 


For  example,  in  the  three-element  array  of  Fig.  11-2 

=  0.0739  +  j0.5621 

12v 

=  -7.075  +  j0.7090  ll:l6 

22v 

=  0.  1605  +  jO.5787 

The  relation  between  the  A  and  B  current  coefficients  is  reduced  to  the 


following  form 
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■•^^2*22v 


4 


21uJ 

u 
Z2u. 

u  > 


j2Nu 
U  : 


where  4',/'^  is  the  element  in  the  second  column  and  i^^  row  of 

u  I  u  ■* 

With  the  simplification  offered  by  11;18,  the  element  currents  are 
recalculated  in  Appendix  XI,  with  the  result 


(1.292  tj6.737)F„^jl 
‘«2<=‘>  =  ><>'\2l-j2.4124Mh^2Mn.ll  -j27.22)F^^^l 
=  ‘°'\2'  (-17.88.jl8.28)F_^^3l 


11:18 


11:19 

11:20 

11:21 


The  validity  of  assuming  that  all  are  negligible  except  *22v  *^7 

comparing  11:19  -  11:21  with  11:9  -  11:11. 
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APPENDIX  VI 

For  convenience,  the  relation  between  the  base  currents  and  the  driving 
voltages  for  ph  =  u/Z  is  repeated  below,  or,  from  6:51 


I*  ] 

u 


2it 


i« 


u 


VI:  1 


The  three-element  array  with  quarter-wavelength  spacing  and  the  main  beam 
in  the  endfire  direction  has  the  following  relationship  for  the  base  currents: 


,(0)' 

1 

zl 

ZZ'"' 

1  =  I,,(0)  . 

-j 

z3'“l 

1 

-1 

VI:2 

» 


The  elements  of  the  matrices  in  VI:  1  are  given  by  6:26  and  6:27  with  6:30 
6:35,  or 


♦l2u  =  -*210  ='»23u  =  =  f 


•*,-3  =-0.4988  +  10.2089,  Bb  =  tr 

13u  31u 


’^'iv  =  =  *3^  =  .  n-  =■  10 

=  *21v  =  *L  =  ■*32v  =  -»•  -  ■i'*-  ■  P»  =  i 


>  VI:  3 


>  VI:4 


^Uv  =^?iv  =  -0-2925  +  j0.1l86  , 


pb  =  Ti' 
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Also 


Zn  -  3 

X  10  .  VI:5 

c^o  dR 

The  inverse  matrix  on  the  right-hand  side  of  VI:  1  is  constructed  from  the 
values  in  VI: 3,  VI:4  and  VI: 5,  with  the  result 


-j(0.4145)x;i0.^ 


(-^,3874-j0.5097) 

(-ai86l-j0.2828) 

(-0.2063+j0.0903) 


(-C.186l-j0.2828) 

(-6.3874-j0.5097); 

(-0.186l-j0.2828) 


-1 


,(-0.206  3 +j  a  090:^1 

(-0.1861 -j  0.2828)) 
(-6.3874-ja5097) 


VI:6 


The  inverse  matrix  VI:6  is  calculated  by  the  same  method  employed  in 
Appendix  I.  First,  all  the  double  products  of  the  individual  elements  of  VI:6 
are  obtained: 


=  40.  5391+j6.  5114 
=  -0.  0458+j0.  1052 

1  w 

^13  ""  0345-jO.  0373 


4  =  1.0446+jl.9O13 

XX  X  w 

♦ll#13  =  1.  3637-j0.47l6 
*1^13  ^  ®^39+j0.  0415 


>  VI:7 


J 


where 


-6.3874-j0.5097,  <*^2=  "O- 186l-j0.2828,  <*^3= -0.2063+ja0903 


.  VI:8 
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The  complete  inverse  matrix  VI:6  is  given  by  VI:7  and  1:3,  with  the  result 


2ir 


L  <ro  dR 


Ztt 


-1 


10^(a3664+jl.5307) 


(a4058+j0.0641) 

(a0098-j0.0186) 

{-aoi4i+jaoo5^ 


(0.0098-ja0186) 

{0.4060+ja0655) 

{a0098-ja0l86) 


(-a0141+j0.005^"^ 
(0.0098- j0.0l86) 
(0.4058+ja064l) 


VI:9 


When  VI:9  is  multiplied  by  the  matrix  I^],  and  the  result  substituted  in  VI:1, 
it  follows  that 


10^(a3664+jL53Cr7) 


(a3459-ja4541)  {-0.1534-j0.3240) 
(-0.16  39 -ja  3196)  (a3928-ja4443) 
.(-a235bfja076l)  (-ai534-ja3240) 


(-0.2356+j0.C76l) 
(-0.16  39 -j  0.3196) 
(0.3459-ja4541) 


VI:10 


The  specified  base  currents  VI:2  are  substituted  in  VI:10,  and  after 
multiplication  the  result  is 


fv  1 
01 

0.2575-j0.3768 

0.6  711 +j  a  2561 

-^02 

>  =  (0.3664+j  1.5  307)xl0‘^< 

-0.444  3-j  a  39  28 

►I  ,(0)  =  10^.' (0)< 
zl  zl 

0.4385-j0.8340  V 

^03 

V.  ^ 

-0.  9055+j0.6836 

-1.3782-jl.l355 

J 

(0.  6711+jO.  2561)1^^(0)1 


10  < 


(0.  8340+j0.  4385)1  JO) 
tZ 

(1.  3782+jl.  1355)y0)^ 


VI:  11 


VI:12 
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The  driving-point  admittances  are  given  directly  by  VI:  12,  or 


=  (13.  01  -  j4.  963)x10'^  mho 
Yo2  =  (9.  393  -  j4.  939)xl0'^  mho 
Yp^  =  (4.  322  -  j3.  560)xl0‘^  mho 


VI:13 


The  quasi  zeroth-order  element  current  is  given  by  6:19  with  the  A 
coefficients  given  by  6:25,  and  the  B  coefficients  given  by  6:49  and  6:50  with 
VI:  13  and  6:25.  With  respect  to  the  individual  driving  voltages,  the  element 
currents  are 


^3.  01  COsPz  -  j 

7,  375cospz  +  2.  4124(sinp|z| 

i  1 

r 

393cos  pz  -  j 

7.351cosPz  +  2.  41 24(sinp  jz|  -  1^ 

1 

•  VI:14 

Z3'“l  =  ‘"'Xs] 

r  1 

322cos  pz  - j  j 

5.972co8pz  +  2.  41 24(sin  pjzj  -  1^ 

} 

• 

The  three-element  currents  are  given  with  respect  to  a  common 
reference  by  expressing  V^^  and  V^^  in  terms  of  V^^  from  VI:11,  with  the 
result 


Vq^  =  (0.0909  +  j0.  7569)Vp^ 

=  (0.3859  -  jl.  8553)Vu2 
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With  VI:  15  in  VI:14,  the  element  currents  with  respect  to  V^^ 


I^^^)  =  10  ^1.826(sin  -  :1)  +6.765cospz  -  j  jQ,2193(sin(3lz|  -  1)  -  9.177cosp^| 

=  10  ^^2^  9.393cos(3z  -  j  ji4124(sinp[z|  -  1)  +  7. 351cosp^^ 

I^(z)  =  10■^^^2|A4759(8in  p|zj  -  1)  -  9.412co8pz  -  j  Ja9309(sinp|z|  -  1)  +  lG.32cosp:jj 


VI:l6 


SRIZ 


VII- 1 


APPENDIX  VII 

The  relation  between  the  specified  base  currents  and  the  driving 
voltages  is  given  by  Z:50  and  repeated  below  for  convenience: 


+ 


1  -  cosph  L  uj  L  u! 


VII:  1 


With  ph  =  3v/4,  the  functions  in  Z:53  -  Z:58  are  evaluated  at 

z  =  0,  except which  is  evaluated  at  z  =  V8.  The  '%i  functions  are  sub¬ 
stituted  in  Z:34  and  Z:35,  with  the  result 


^llu=  ^ZZu=  = 

^IZu  ^  *Zlu  =  ‘'Z3u  =  ‘*’3Zu  ^  0.9148431.4058  .  pb  =  f 

^  =  0.9861 -jO.  5075  .  pb  =  n 

1 3u  31u 

=  0.  3Z15-jl.  8139  ,  Sl.=  10 

llv  ZZv  33v 

=  -0.647Z-j0.9839  ,  pb=  ^ 

IZv  Zlv  Z3v  3Zv  ^  Z 

,  =*^^1  =  -0.  7140+j0.  36Z5  ,  pb  =  T7 

13v  31v 


j'VIPZ 


VII:  3 


Also, 


Ztt 


-J3.  3Z19xl0‘^  ,  -j  — =  j3oi.  034  VII:4 

2tt  • 


The  inverse  matrix  on  the  right-hand  side  of  VII:  1  is  formed  with 


VII:Z,  VII:3  and  VII:4,  or 


SR12 


VII- 2 


sinph  r^-i 
1  -cosph  L  uJ 


-1.9202-ja7427 

-0.2683-j0.40l6 

-0.3056+j0.1523 


-a2683-j0.40l6 

-l.9202-j0.7427 

-a2683-j0.40l6 


-0.3056+j0.1523 

-0.2683-j0.40l6 

-1.9202-ja7427 


-1 


VII;  4 


The  inverse  matrix  VII:4  is  found  by  the  same  method  employed  in  Appendix  I. 
First,  all  the  double  products  of  the  individual  elements  of  VII:4  are  obtained: 


=  3.  1356+j2.  8522 
=  -0.  089  3+j0.  2154 

1  u 

=  0.  0702-j0.  0931 


^ll‘**12  ""  0.  2l69+j0.  9705 
4ii4i3  =  0.6999-j0.0654 

*12‘*'l3  ^ 


VII;  5 


where 


4^^  =  -l.9209-j0.7427;  4^^  =  -0.  2683-jO.  4016  ;  4^  ^  =  -0.  3056+j0.  1 523 

_  VII:6 


The  final  form  of  the  inverse  matrix  VII:4  is  given  by  VII:5  and  1:3,  with  the 
result 


1  -  cos  ph  LuJ 


{-0.05923+j0. 10019) 


-2.3725-j0.6l72 

0,0836+jl.0783 

0.9245-j0.35l7 


0.0821+jl.07l5 
-2.  l843-ja9847 
0.0821+jl.0715 


0.9245-j0.35l7 
0.0836+jl.0783 
-2.3725 -j0.6l72 


VII:  7 


VII- 3 


SR12 

When  VII:7  is  multiplied  by  the  matrix  and  the  result  substituted  in 

VII:  1  it  follows  that 


^1=  10^(-0.1767-j0.1044) 


-2.3725-j0.6712 
0.0836+jl.0783 
0. 9245 -j  0.3517 


0.0821+jl.0715 
-2. 1843-a9847 
0.0821+jl.0715 


0.9245 -j  0.3517 
0.0836+jl,0783 
-2. 3725 -j  0.617  2 


VII  :8 


where  for  the  endfire  array 


.  VII:  9 


With  the  base  currents  specified  by  V£I:9.  the  driving  voltages  are  given  by 
VII:8,  or 


r 

V 

-2.2255 -jO.  3476 

0.3570+j0.2937 

01 

^  ^02 

>  =  loio.]%7.j0.10r44)I^^(0)  < 

-0.9847+jl.9257 

>  < 

0.4020-j0.3743  V 

V 

4.3685+j0.1834 

-0.5408-j0.4885 

l;o3j 

^  ’J 

^1 

(0.  3570+j0.  2937)1^^, (0) 

o 

> 

>  =  10^< 

(0.  3743+jO.  4020)1^0)) 

(0.5408+j0.  4885)1  {O) 

The  driving-point  impedances  and  admittances  are  given  by  VIIill,  thus 


SRI?. 


VII -4 


=  357.  0+j293.  7  ohms 

=  (1.6706-jl.  3744)xl0'^  mho  > 

Z  =  374.  3+j402.  0  ohms 

U  u 

=  (1.  3324-jl.  2402)xl0'^  mho 

Z^3  =  540.  8+j488.  5  ohms 

Y^3  =  (1. 0183-j0.  9198)x10'^  mho 

The  quasi  zeroth-order  current  is  given  by  2:27,  2:45  and  2:47  with  2:40 
and  VII:11.  Hence,  with  respect  to  the  individual  driving  voltages,  the  currents 
are  given  by 


I  fo.  9786F  -j{3.3219M  -  0.  5708F  )lxl0 

zl  01  •-  oz  oz  oz  J 

I  To.  7805F  -j{3.  3219M  -  0.6495F  )lxl0' 

z2  02  L.  oz  oz  oz  J 

I  =  V.,  fo.  5965F  -j{3.  3219M  -  0.  8371F  )lxl0‘ 

z3  03  i-  oz  oz  oz  J 

where 

^oz  ""  T  ■ 

F  =  cos  flz  +  0.  7071 
oz 

The  three-element  currents  are  given  with  respect  to  a  common  refer¬ 
ence  by  expressing  and  in  terms  of  from  VII:  10,  with  the  result 

V^^  =  (0.  1115+j0.834l)V^2 

^02  =  ^02 

V^3M-0.I148-jl.32l6)VQ^ 


>  VII:13 


J 


^  VII:  14 


>VII:15 


SRIZ 


VII- 5 


With  VII:15  in  VII:13,  the  element  currents  with  respect  to  are 

I  ,(z)  =  V  ^  r?.  7618M  .  0.  3669F  -  j{0.  3704M  -  0.8798F 

zl  02  L  oz  oz  oz  oz  J 

I  o(z)  =  r  0.  7805F  -j(3.  3219M  -  0.6495F  )lxl0‘' 

z2  02  (_  oz  oz  oz  J 

I  o(z)  =  V„,r-4.  3902M  +1.0378F  +j(0.3818M  -  0.  8844F  )!  xl0‘ 

z3  02L  oz  oz  oz  ozJ 


VII:  16 


SR12 


VIII- 1 


APPENDIX  VIII 

The  general  relation  between  the  base  currents  and  the  driving  voltages 
for  the  N-element  array  is  given  by  2:50 


1^,(0) 


1 1  -  cos(3h 


=  J 


Zir 


/o 


sin  (3h 


1  -  cos^ 


VIII:  1 


where  for  the  parasitic  two-element  couplet 


VIII:  2 


With  ph  =  3tt/4  and  pbj2  =  the  values  of  the  matrix  elements  are  given 

in  Appendix  VII,  or 


f#l  = 

-5.  4124+j2.  5861 

0.  9148+jl  .4058” 

L  uJ 

0.  +148+jl.4058 

-5.  4124+j2.  586 1_ 

VIII:  3 


_ 1, 

1  -  cosph  u-l 


-1.9202-j0.  7427 
-0.  2683-jO.  4016 


-0.  2683-j0.40l6 
-1. 9202-j0.  7427 


VIII:4 


The  substitution  of  VIII:3,  VIII:4  and  VIII:2  in  VIII:1  yields  the  following 
result: 


-j3.3219xl0~^Voi 
23.74,61 -j30,56<'-i) 


-5.4124+j  2.5861 
-0. 9148- j  1.4058 


-0.9148-j  1.4058' 
-5.4124+j  2. 5  861 


'-1.9202-j0.7427l 

-0.2683-j0.40l6 


VIII  :5 
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vin-2 


where 


ffo  o  dR 


After  performing  the  matrix  multiplication  in  VIII:?,  the  result  is  given  by 


K(0) 


1.7071^ 


ri2.  0044-j0.  2014 
J  =  lO'^V  (6.  7773-j5.265lW 
^  [_3.  2082+j4.  8585 


VIII;  7 


=  10* 
11.70711  01 


[0.80299-j0.  64778 
JD.47323+jO.  16036J 


VIII;  8 


The  two  base  currents  I  (0)  and  I  (0)  are  given  directly  by  VIII;8.  Thus, 

Z  X  Z  w 


=  10  V 


01 


1. 3708-jl.  1058 


0.8079+j0.  2738J 


VIII;9 


The  B  coefficients  of  the  element  currents  are  given  by 


B. 

1 


1 


1  -  cosfJh 


-  jA.sinfh^ 


VIII;  10 


where 


-jA^ sin ph  =  j2.  3489  VIII;11 

jA^  =  -j3.  3219  X  10 


VIII;  12 


SR12 


VIII;  3 


and 


jA^  =  0 

The  element  currents  are  determined  from 


I  , 
zi 


(z) 


jA,M  +  B.F 
1  OZ  1  oz 


VIII:  13 


VIII:  14 


where  the  A  coefficients  are  given  by  VIII:  12  and  VIII:  13,  and  the  B  coef¬ 
ficients  are  given  by  VIII:  10  with  VIII:9,  whence 


I  ,{z)  =  10"^V-,  r0.8030F  -j(3.3219M  -  0.  7282F  )1 

zl  OIL  OZ  oz  oz  -J 

I  Jz)  =  lO' V,  (0.4732F  +j0.l6O4F  ) 

z2  01'-  oz  oz- 


VIII:  15 


The  rigorous  solution  for  the  same  couplet  described  above  follows 
from  the  symmetrical  component  solution,  or 


(z)  = 


m 


(M  +  'r'''‘(h)F  ) 

oz  oz 


VIII:  16 


dR  o 

where  for  the  couplet,  m  =  0  corresponds  to  driving  both  elements  in  phase, 
and  m  =  1  corresponds  to  the  phase  reversal  case.  The  T'^'(h)  functions  are 
given  by 


T^(h)  =  ^12v 

'^’llu  '*’l2u 


VIII:  17 


where  the  upper  sign  corresponds  to  m  =  0,  the  lower  to  m  =  1.  The  numerical 


SR12  VIII-4 

values  of  the  1^))  functions  are  easily  calculated  from  the  #  functions  given 
in  Appendix  VII,  or 


T®(h)  =  -0.  2367  +  jO.4646 
T^h)  =  -0.  1720  +  jO.  0994 


.  VIII:  18 


The  phase  sequence  currents  are  given  by  VIII:  16  with  VIII:  18; 


I®  (z)  =  I^,(z)  =  lO'^V®  fl.5434F  -j(3.3219M  -  0.7863F  )] 

z2  zl  >-  oz  oz  oz  J 

VIII:  19 

-I  (z)  =  l\(z)  =  lO'^V^fo.  3301F  -j(3.  3219M  -  0.  5714F  )! 

z2  zl  L  oz  oz  oz-l 


where 


0  _  1  _  I 

~  2  01 


VIII:  20 


The  phase  sequence  currents  are  combined  to  give  the  element  currents  for 
the  couplet  with  V^^  =  0.  thus 


I  (z)=|r  +4l  =  lO'^V-,  ro.9368F  -j(3.3219M  -  0. 6789F  )] 

zl  2  z  2  z  Oil-  oz  oz  ozJ 

VIIL21 

1  ,0  11 


I  -(z)  =  -  I  -71=  10"  V-,  r0.6067F 
z2  2  z  2  z  01  L  oz 


+  jO.  1075F^  j] 
oz  J 


SR12 


IX- 1 


APPENDIX  IX 

The  physical  dimensions  of  the  couplet  are 


pb  =  1  .  10  =  In(^)  =  ln(-^) 

■^12  aj  &2 


IX:  1 


ph^  =  1.4  ,  ph 


2  2 


and  the  driving  voltages  are 


^01  =  °  •  ^02 


IX;2 


The  forms  of  the  'P  functions  which  are  needed  to  compute  the  ♦  matrices 
are  given  below 


IX  :4 


12du  1  -  cosph^ 


[c^(h^-,0)  -  C^di^.hj)] 


IX:5 


W-i' 


21du 


IX  :6 


=  -  co.phj  IX:7 


IX:8 


=  C  (h_,0)  -  C  (h„.h,) 
22du  a  2  a  2  2 


IX  :9 


SR12 


IX-2 


izdv  =  -  K‘V®-\'VV]] 


■^12v'N'=W‘'i>-  W"!* 


22dl 


=  Im  -  fa'V‘«-\'VV]) 


■^22v<“'  =  \<''2’  O* 


^  0.7071Ra-i-,^^„(0)^  0.2929Re-y,^^„(;) 
-0.4142 


12du 


■Re:il2du(0)  -  Re-tiWf) 


-0.4142 


•4dR  =  -«'{[\"'2’<»-®a<''2->'2']  -  K'>'2’'»-=al»2->'2>]} 


The  following  sine  and  cosine  integral  function  values  are  needed 
evaluate  the  preceding  >4^  functions  (N.B.  £X.=  10  for  both  elements.): 


E(h  ,0)  =  9.  1108-j2.5124 

3i  X 

E^(h  .  0)  =  0.  7458-j2.C968 
b  1 

E^(h^,h^)  =  4.  1868-jl.  8321 

^b-^1’^2^  =  -0.  1701-jl.  3249 

E^(h^  0)  =  8.  8836-j2.  7784 

Eb(h^,  0)  =  0.6950-j2.  3126 

E^(l^,l^  =  4.  1786-jl.  8514 

Ebih^.hi)  =  0.  04l3-jl.6607 

IX:  10 

IX;  11 

IX:12 

IX;13 

IX:  14 

IX:  15 

IX;16 

IX:  17 

to 


SRl^ 


IX- 3 


C  (h  ,0)  =  8.6011-jl.8321 

ct  1 

""  1 -4644-31.3205 

3495-jl.  8514 

^a^^’^2)  "  0-5737-jl.  1976 

S  (h^,  0)  =  1, 8514-11.6851 
a  2 

S  (h  ,h  )  =  3.7364-31.1549 


Cj^(h^,  0)  =  0,  6670-31. 5335 
C|^(h^,h^)  =  -0.1809-30.9851 

^b^^2’  5490 

CJh^.h^)  =  -0.  0431-31.  1197 

S^(h^,  0)  =  0. 2278-31.  3956 
b  2 

Sb(I^,hi)  =  0.  0951-31. 0247 


IX:  18 


The  components  of  the  matrices  in  9:35  and  9:36  are  calculated 
from  9:37  -  9:40  with  IX:2  -  IX:18,  {N.  B.  Since  =  0,  it  is  not  necessary 
to  calculate  every  member  of  all  the  matrices.  )  with  the  result 


‘^1.  =  '*','^.  =o;2179 
21c  1 2du 


*llu  =  =  0-”71+j0,9580 


<iu  =  <  ^Idu  "  i‘"’'*21du>  -  '^21u<'"  = 


IX  :1  ? 


Klu  -  W*'lh.  ’^2u<S'  =  »■  >925+jl.  0764 


♦•>5  =  %->j  '  (“I  =  “■  l605+i0.  5784 

7.2u  22du  22u 


:*7“7v  =  ,1i2v'“lr  ’O’" 


From  9:33  note  th?t  with  =  0,  iA^  =  0  ?nd  the  vplue  of  -3A^  is  given  by 


SR12 


IX:4 


-jA 


2 


Bl,tl2du 

'*22dR 


When  the  value  of  -jA  from  IX; 20  is  substituted  into  9:36,  the  E 

Ct 


are 


given  by  the  4olution  of 


ray = [<1 


.i _ ^2 _ V  -  ^ 2dn 


22dR 


which  may  be  rearranged  in  the  form 


aU  +4^  lUdu 

^u  ” 

12u 

LX  -  ; 

fU  tUdu. 

^2'^'^22dR 

^u 

22u 

J  <5o  22dR 

12v 


‘'02  \ 


22v, 


When  IX:22  is  solved  for  the  B  coefficients,  the  result  is 


2tr 


■J 


■V, 


22dR 


02 


^l|du 

♦rzu 

-1 

r  N 

12v 

< 

tigdu 

21u  22vf'  n 

L  22dR 

22v 

The  components  of  the  matrices  in  IX: 2 3  are  given  by  IX:  19.  or 


IX:  20 


coefficients 


IX:  21 


IX;  22 


IX:23 


j2.4124xlO'^V, 


02 


0.5422+j0.9762  0.1925+j  1.0764 

0.1263+j0.9519  -0.6880+jl.2l87 


r\  -1 


[0.l605+j0.5784 

-7.0754+j0.7D90j 


The  inverse  matrix  on  the  right-hand  side  of  IX: 23  is  given  by 


*U  .u 

l5lu  22v^ 


=  (-1.32  31 +j  0.7  760) 


-0.6880+jl.2187  -0, 19  25 -j  1.0  764 


-0.1263-j0.9519  0. 54  22 +j  0.9  76  2 


With  IX;25  in  IX:24,  the  B  coefficients  are  given  by 


B 


B. 


-20.  78+jl7.  81[ 
Il4. 06-j25. 39. 


xlO'^V 


02 


The  coefficient  of  the  M  current  on  element  two  is  given  by  IX:20  with 

oz 


or 


-JA.,  =  r-j2.  4124  -(20.  78  +  jl7.  81)(0.  03154)1x10' , 

c  L  ^  [jC 


-JA  =  (-0.6554  -  j2.4686)xl0'  V 


^  IX:24 


IX:  25 


IX:  26 


IX:  26, 


IX:  27 


The  current  on  each  element  is  given  by  9:29  and  9:30  with  IX:26  and  IX:27: 


SR12 


IX-6 


I  ^(z)  =  [(-0.6554  -  j2.468)Mh^2  +  (14.06  -  j25.  39)F^^2]xl  O"  ^ 


(jRlZ 


X-1 


APPENDIX  X 

The  physical  dimensions  of  the  couplet  are 


Pb 


12 


=  1 


XI  =  10  =  In  ^ 

^1 


ph^  =  1.8  ,  =  I 


X:1 


The  driving  voltages  are 


V 


01 


-  0  , 


X:2 


The  general  forms  of  the  functions  are  given  in  Appendix  IX.  The  additional 
values  of  the  sine  and  cosine  integral  fiinctions  required  for  this  couplet  are 
listed  below. 


C  (h,,h,)  =  -1.4153-jl.027V 

ail 

C  (h.  .0)  =  8.  1053-jl.  8219 
a  1 

Cj^(h^,h^)  =  -0,  3508-j0.  5527 
Cj^(h^,0)  =  0.  6742-jl.  5254 

Cb(h2,hi)  =  -0.  2983-j0.8329 
S  (h  ,  h  )  =  -0.  0840-j0.  7855 


E  (h,,h,)  =  3.6910-jl.  8219 
all 

E^(h^,  0)  =  8.  6050-j3.  0116 

^b^^l’^2^  =  -0.  0448-jl.6696 
Ej^(h^,  0)  =  0.6l58-j2.4992 

^b^^2’^l^  =  -0.  2784-jl.  2670 
S|^(h^,0)  =  -0.  1496-jl.  5804 


X:3 


Since  V 


01 


0  and  h^  >  h^, 


the  sinusoidal  current  on  element  one  is 


zero,  and 


SR12 

jA^  =  0 


X-2 


-jA 


2 


X:4 


The  elements  of  the  4  matrices  for  hj^  h^  (i.e.  9:41  -  9:44)  j  are  given 

below: 


u  =  =  ->•  3920+jl.6339 


=  *9.^  -  *9.  (0)  =  0.  36lO+jO.  9320 

21u  21du  21u  ■' 


**22t  ""  ^22dR  "" 

*12v  ”  '^IZv  ■  ^o^^l^^'^Udv  •^^"*'^12dv^  "  5621 


X:5 


4v  =  -  ^'^ZZdl  -  -7.0754+30.  7090 


With  X:4  in  9:36,  it  follows  that 


-J 


<? 


2it 

vi.h 

o  22dR 


X;6 


The  inverse  matrix  on  the  right-hand  side  of  X:6  is  computed  from 


SR12  X-3 

the  values  given  in  X:5,  with  the  result 


(-0.  0176+j0.  2993) 


-0.  6880+jl.  2187 
-0,  36l0-j0.  9320 


-0.  1722-jO.  9655 
-1.  3920+jl.  6339 


X:7 


The  inverse  matrix  X:7  is  substituted  in  X:6,  and  after  matrix  multiplication 
the  B  coefficients  are  given  by 


(a7220+j0.0425)y 


02 


1.  I67l+j6.4126 
9.  1877-jl2.  8192J 


xlC 


-3 


ro.5701+j4.6794' 
1.7.  1783-j8.  8649., 


X:8 


The  element  currents  are  given  by  9:29  and  9:30  with  X:4  and  X:9i  or 

I  ,(z)  =  10‘ V,(0.5701+j4.6794)F  , 

Z 1  \JC  oz  1 

I  _(z)  =  10" V  r-j2.4124M^  ,  +  (7.  178-j8.865)F  .I 
z2  02  L  •'  oz2  oz2J 


X:10 


SRI  2 


XI- 1 


APPENDIX  X.I 

The  relationship  between  the  B  and  A  coefficients  for  the  three-element 
Yagi  array  of  Figure  11-2  is  given  by 


(b)  =  Kf  m  {iA2 


where  from  11:7  and  11:8,  for  =  1.8,  ph^  =  Tr/2,  ph2  =  1.4,  and 
pbi2  =  1  =  pb^^ 


=  -1.  3920+j  1.6339 

1  lu 


=  0.  36l0+j0.  9320 

21u 


^r2u  ""  ^ 


^22u  ""  -0.688C+jl.  2187 


^13u  ■  ^•'^326+jO.  3073 


4“  =  0.1323+30.9295 

23u 


=  0.8109+30.  5728 
31u  . 


*12v  ^  0'^^9+30.  5621 


=  0.1925+31.0764 

32u 


=  -7.  0754+30.  7090 


=  0.5371+30.  9580 

33u 


=  0.  1605+30.  5784 


The  inverse  of  the  matrix  is  symbolically  given  by 


SR  12 


XI-2 


where 


-*21*33+ *23*31 

=*21*32  -*22*31  = 

=  -^^2^33  ^^13^32  "  0-5941-jO.  lObV 

*''=*11*33  -*13*31 

*"=  -*11*23+^13*21  = 

*‘'  =  *12*23  -*13*22 

4^^  =  ^‘^3l‘*’l2  ""  1-56V5-Ij2.  0344 

*“=*11*22  -*12*21  =  -“•'^*‘-■'^•^■'’2 
A=  ■'■‘^21^^^  ^'*'31’*^^  ^  1.5035-j0.  1332 


XI:4 


With  XI:4  and  XI: 2  in  XI:  1,  the  B  coefficients  are  related  to  -jA^  by 


B. 


'2  = 


-0.5620-j0.3258 

0.2720+j0.0372 

-0.7935-j0.8l43 


0.594l-jC.1069 

-2.3836-jl.0047 

1.5737-l-jl.5918 


-0.2  0  97-j0.0782 
1.5695+j2.0344 

-0. 1646 -j  3. 27 92  I 


0.0739-l-j0.5621 
<'-7.0754+j0.7090)>  XI5 
ai605-(-j0.5784 


The  result  of  the  matrix  multiplication  in  XI: 5  is 


SR-12 


XI- 3 


'-2.  6704+j0.  3049" 

®2 

10.5901+j5.4671  1 

^-5.  9350-j8.  0246_^ 

XI:6 


whe  re 


-jA^  =  -j2.4124V^2*10 


XI;  7 


With  XI:7  in  XI:6  it  follows  that 


-3 

^  “V 

0.  7355+j6.4421 

1^ 

o 

=*s3 

O 

> 

II 

13.  l888-j25.  55  1 

-19.  3585+jl4.  0764; 

XI:8 


From  XI;8  it  follows  that  the  element  currents  are  given  by 


(0.  7355  +  j6.4421)F  ,1 

ozlj 


'z3W  =  “>'\2r 


(-19. 36  +  jl4.  08)F^^2  J 


When  it  is  assumed  that  is  large  compared  with  and  , 

22v  or  j2v  32v 


then 


SR12 


XI-4 


- 

r 

■^l' 

0. 5941-jO.  1069  -- 

... 

li 

1 

> 

1 

-2.  3836-jl.  0047  -- 

< 

-7, 0754+j0.  7090  > 

B3 

1.5737+31.5918  -- 

or 


'«i^ 

1. 2921+j6.  7374I 

-3 

I 

>  = 

^02- 

1  1.  1067-327.  22 

.^3. 

-17.  88+3I8.  28 

1  y 

The  element  currents  for  this  case  follow  from  Xl:ll,  or 

=  ‘“■\2  r  (1.2921  +j6.737)F^^J 

‘z3'''  =  “>'\2  [  (-17.881.  jl8.28)F^^3] 


Xl:10 


XI:11 


XI:12 
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